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ABSTRACT 

A  digital  computer  program  is  presented  capable 
of  determining  the  in- and/or  out-of-plane  vibration 
frequencies  of  a  single  plane  piping  system0   The 
analysis  is  made  utilizing  the  method  of  transfer  matrices, 
a  discussion  of  which  is  included »   The  program  accepts 
branched  systems  with  or  without  intermediate  supports 0 
The  effects  of  both  shear  deflection  and  rotational 
inertia  are  considered  although  either  or  both  may  be 
neglected  at  the  user*s  discretion.,   Distributed  mass  is 
generally  treated  as  such,  although  provision  is  made 
for  alternate  treatment  by  lumping  mass „ 

Several  typical  oiping  systems  are  analyzed  using 
these  various  program  options,  and  the  -results  obtained 
are  compared. 
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CHAPTER  I 

INTRODUCTION 
Id  General  Remarks 

The  dynamic  behavior  of  piping  systems,  particularly 
their  response  to  periodic  and/or  transient  excitation, 
is  a  subject  of  considerable  interest  and  growing 
pertinence  in  the  design  of  such  systems 0   The  reasons 
for  this  are  readily  discernable  when  the  unprecedented 
reliability  demanded  of  nuclear  power  plant  piping,  or 
the  consequence  of  a  piping  system  failure  in  an  orbiting 
space  vehicle,  are  consideredo   An  increasing  number  of 
piping  systems  are  being  designed,  of  which,  like  those 
mentioned  above,  a  long  service  life  is  demanded,  and 
to  which  access  is  limited  or  completely  denied  once  the 
system  Is  placed  in  service..   Clearly,  no  factor  which 
might  contribute  to  the  failure  of  such  a  system  can  be 
ignored.   Fatigue  brought  about  by  periodic  excitation 
at  or  near  a  system  resonance  is  one  such  factor „   Cata- 
strophic failure  caused  by  a  periodic  disturbance  of 
short  duration  and  high  amplitude,  such  as  an  earthquake, 
is  another,* 

With  regard  to  the  latter  case,  It  has  come  to  the 
writer's  attention  since  beginning  this  work,  that 
Japanese  engineers  are  giving  considerable  and  close 


attention  to  this  particular  problem*  Their  first  step 
in  approaching  it  is  to  determine  the  natural  frequencies 
of  the  system  of  interest 0 

Fortunately,  or  perhaps  fostered  by  necessity, 
development  of  the  "tools"  needed  to  carry  out  accurate, 
reliable,  and  practicable  vibration  analyses  of  complex 
systems,  namely  analysis  by  matrix  methods  with  the  aid 
of  high  speed  digital  computers,  has  kept  pace  with  the 
increasing  demand  for  their  employment „ 

1.2  Scope  of  Work  Presented 

A  method  for  determining,  analytically,  the  in- and 
out~of-plane  undamped  natural  frequencies  of  planar 
piping  systems  has  been  developed  and  is  presented 
herein,.   It  employs  the  method  of  "Transfer  Matrices"1 
to  construct  an  accurate  mathematical  system  model  and 
determine  system  natural  frequencies,  and  a  digital 
computer  to  carry  out  the  inherent  voluminous  numerical 
computations o  VIPIPE,  the  digital  computer  program  pre- 
sented in  appendices  A  through  C  constitutes  the  embodi- 
ment of  the  method „ 


*The  theory  and  method  of  Transfer  Matrices  is  dis- 
cussed in  sufficient  depth  for  the  purposes  of  this 
report  in  chapter  2„   A  more  comprehensive  and  detailed 
treatment  of  the  method  may  be  found  in  refo  L 


This  method  is  a  first  but  important  step  toward 
the  ultimate  objective  of  analysis  of  three  dimensional 
systems  of  unrestricted  topology 0 

Attention  has  been  confined  to  the  two  dimensional 
or  planar  case  because  the  theory  underlying  the  treat- 
ment of  three  dimensional  systems  using  Transfer  Matrices 
is  not  presently  available  in  the  literature,  and  because 
it  is  felt  that  the  need  for,  and  utility  of,  a  treat- 
ment of  the  two  dimensional  case  fully  warrants  the 
development o 

Although  in  its  present  state  of  development,  the 
theory  of  Transfer  Matrices  permits  the  consideration 
of  any  planar  piping  configuration,  certain  restrictions 
are  imposed  on  the  configurations  which  may  be  treated 
by  the  method  herein  presented, by  the  limitations  of  the 
program  VIPIPEo   These  limitations  are  discussed  in  detail 
in  section  103" 

The  philosophy  of  VIPIPE  has  sufficient  generality 
to  make  the  program  readily  adaptable  to  consideration 
of  the  three  dimensional  case  as  soon  as  the  underlying 
theory  is  developed • 

1  <>3  Assumptions  and  Limitations 

The  mathematical  models  used  to  describe  the 
elastic  properties  of  a  system  and  its  components  are 
faithful  within  the  framework  of  the  following  con- 
siderations : 


au   The  material  of  which  each  component  is  made 
is  considered  to  be  homogeneous,  isotropic,  and 
linearly  elastic c 

bo   The  system  is  conservative  (no  damping) „ 
c0   Distributed  mass  is  treated  as  such  except  in 
curved  sections  for  which  a  lumped  mass  model  is 
usedo   For  curved  sections,  sufficient  sub- 
sectioning  is  carried  out  so  as  to  reduce  error 
to  an  economically^acceptable  level » 
do   Spring  hangers  are  considered  to  consist  of 
from  one  to  three  mutually  perpendicular  linear 
springs,  and  from  one  to  three  torsional  springs 
which  act  about  mutually  perpendicular  axes 0   When 
a  spring  hanger  consists  of  both  linear  and  torsional 
springs,  the  linear  springs  are  normally  considered 
to  act  along  the  same  mutually  perpendicular  axes 
that  the  torsional  spring  act  about 0   A  more  gen- 
eral case  may  be  treated  with  relative  ease  and 
is  a  subject  of  discussion  in  appendix  B,  section  Bo96° 
Although  the  theory  would  permit  determining  any 
number  of  natural  frequencies  for  any  planar  system 
having  any  number  of  components,  branches,  and  hangers, 
program  VIPIPE,  described  herein,  is  limited  to  the 
treatment  of  systems  having  a  maximum  of  fifty  components, 


^Appendix  B,  section  Bo93. 


twenty  branches,  and  twenty  hangers 0   Further,  only 
single  filament  branches  from  a  main  member  may  be 
treated,  i.eo,  all  branches  must  emanate  from  the  same 
continuous  member,  and  a  branch  may  not  emanate  from 
a  branch o 
l*k     Notation 

[_  ~]  matrix 

x,y,z    cartesian  right-handed  coordinate  system 

u,v,w    displacements  in  the  x,y,z  directions 
respectively 

S.^lx^   rotations  about  the  x,y,z  axes  respectively 


d 

displacement  vector 

£ 

Young 9s  modulus 

£J 

bending  stiffness 

G 

shear  modulus 

GJt 

torsional  stiffness 

G 

coordinate  transformation  matrix 

I 

unit  matrix 

Jt 

polar  second  moment  of  area 

L 

length 

M 

bending  moment 

N 

normal  force 

0 

null  matrix 

P 

internal  force  vector 

H^2 

square  submatrlces 

S 

spring  matrix 

T  torque 

U  transfer  matrix 

Uj*  transfer  matrix  element 

V  shear  force 

v  unknown  vector 

Z  state  vector 

L,R  subscripts,  left  and  right  side  of  com- 
ponents respectively 


CHAPTER  II 

TRANSFER  MATRICES 

2„1   The  State  Vector  and  Transfer  Matrix 
a.   State  Vector     The  state  vector  at  a  point  i  of 
an  elastic  system  is  a  column  matrix  having  sufficient 
elements  to  completely  describe  the  instantaneous  dis- 
placement and/or  change  in  orientation  of  that  point 
from  its  quiescent  position,  as  well  as  the  force  and 
moment  system  existing  in  the  member  at  the  same  point 
at  the  same  instant,  i<,e0,  those  forces  and  moments 
which,  if  a  cutting  plane  were  passed  through  the  point 
at  the  instant  of  interest,  would  have  to  be  applied  to 
each  cut  face  to  prevent  relative  motion  between  them0 

For  perfectly  general  structural  systems  the 
notion  of  a  state  vector  is  a  rather  sophisticated 
concept.   In  any  properly  constituted  state  vector 
however  the  number  of  elements  pertaining  to  forces  and 
moments  is  equal  to  the  number  pertaining  to  displacement , 

In  the  case  of  a  planar  piping  system,  the  state 
vector  has  six  elements,  three  pertaining  to  displace- 
ment, and  three  pertaining  to  force 0   The  state  vectors 
for  the  in- and  out-of-plane  cases^respectively ,  as  used 


^-Henceforth  the  in-plane  case  will  be  referred  to 
as  Case  1,  the  out-of-plane  case  as  Case  20 


in  VIPIPE,  were  chosen  to  be 


-1P= 


^oop- 


(2.1-1) 


The  elements  of  these  state  vectors  are  the  magnitudes 
of  vector  quantities  of  the  same  name  which  derive  their 
orientations  from  those  of  the  x,y,z  axes  of  a  local 
right-handed  cartesian  coordinate  system  (figo  201-1), 
associated  with  the  particular  point  of  interest . 


centre  Id  a  1    axis 


Fig0  2ol-l  Coordinate  system  and  vectorial  representa- 
tion of  state  vector  elements „   (Interpret  double  shaft 
arrows  by  the  right-hand  rule0) 


These  local  coordinate  systems  are  constructed 
such  that  the  plane  of  the  quiescent  system  and  the 
xy  plane  of  the  coordinate  system  coincide 0   The  local 
x  axis  is  coincident  with  the  local  tangent  to  the 
quiescent  system  centroidal  axis.   No  reference  is  made 
to  a  common  fixed  coordinate  system 0   In  establishing  the 
coordinate  system  for  the  purpose  of  employing  VIPIPE, 
the  longest  continuous  run  of  pipe*  or  the  one  from 
which  all  branches  emanate  if  they  are  not  one  and  the 
same,  is  designated  as  the  main  member,  one  end  of  which 
is  chosen  as  the  starting  point „   This  end  is  called 
the  left  end  and  Its  state  vector  is  symbolized  ZL  0 
The  positive  x  direction  is  taken  as  away  from  the  "left" 
end  and  toward  the  "right"  end  along  the  centroidal 
axis  of  the  main  member.,   The  positive  z  direction  is 
perpendicular  to  the  quiescent  plane  of  the  system  and 
invariant,  and  this  direction  is  called  downward  <,   The 
y  axis  is  chosen  so  as  to  complete  a  right-handed  car- 
tesian coordinate  system*, 

be   Transfer  Matrix     Let  us  consider  the  relationship 
between  the  state  vectors  at  the  left  and  right  ends  of 
a  massless  straight  section  of  pipe0   We  begin  by  con- 
sidering the  situation  occurring  in  Case  2  in  which 
both  bending  and  twisting  (about  the  centroidal  axis) 
occur,   The  relationships  between  the  left  and  right 


end  displacements  are  given  by? 


V 


TL 


+  \i+\w7+\-kr 


L 


The  force  relationships  are  given  bys 


VL  2EJ 


tl 


ML  +  VLL 


Vr  = 


(2„l-2) 
(2d-3) 
(2ol^) 

(2ol-5) 

(2.1-6) 

(2d-?) 


or,  in  matrix  notation 


T 

wy 


L 

GJt 
1 

0 

0 

0 

0 


0 

0 

"S 

0 

0 

T 

L2 

2EJ 

l3 
6EJ 

— w 

L 

1? 

zej 

"¥ 

i 

L 

M7 

0 

i 

Vz 

The  six-by-six  matrix  relating  the  state  vectors  of 
the  left  and  right  end  of  the  pipe  is  termed  a  transfer 
matrix o   It  is  the  transfer  matrix  of  the  elastomechical 
system  constituted  by  the  pipe  and  its  boundary  condi- 
tions o   Similar  six-by-six  transfer  matrices  may  be 
developed  for  point  masses,  curved  sections  of  pipe, 


10 


and  elbows,  etc0   Many  of  these  are  contained  In  the 

1 

transfer  matrix  catalog  Included  in  refo  1,  and  are 

used  in  VlPIPho 

Since  the  application  of  interest  is  to  the  deter- 
mination of  the  natural  frequencies  of  piping  systems, 
it  would  be  instructive  at  this  point  to  illustrate  the 
origin  and  nature  of  frequency  dependence  In  transfer 
matrices 0  To  that  end,  let  us  construct  the  transfer 
matrix  of  a  point  mass,  again  considering  Case  2C   To 
begin,  note  that  throughout  the  application  considered 
herein  we  are  concerned  with  the  isochronous  vibration 
of  a  system  at  a  single  frequency CO.   The  tacit  assump- 
tion has  been  made  that  all  quantities  varying  with 
time,  l.e0,  all  of  the  elements  of  the  state  vectors, 
do  so  with  the  form;  constant*  coswt,  and  have  exhibited 
only  the  constant  or  amplitude,  while  suppressing  the 
variable  factor • 

The  transfer  matrix  of  a  point  mass  may  be  constructed 
by  noting  that  deflection,  twist,  slope,  and  moment  are 
constant  and  continuous  across  the  zero  length  section 
which  a  point  mass  constitutes «   The  vibrating  mass  does 
introduce  a  discontinuity  in  shear  however,  which  may  be 
determined  by  considering  the  free  body  diagram  of 


1See  App0  F. 
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Fig,    2,1-2,    and   the   following  relationships; 
w(t)    ■  w»coso>t 

,2  2 

d  w(t )    =  "WcocostDt 

dt2 


(2,1-9) 


From  F=ma: 


V3RoosGJt-VsLcoscat=m<*  w(t) 
dt 


(2,1-10) 


and  substituting  (2,1-9)  into  (2,1=10)  we  have; 

VzH  =  VzL  +  mo2(-w)  (2,1-11) 

and  are  ready  to  construct  the  transfer  matrix  Up  of 
a  point  mass.  Hence 


1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

_° 

0 

2 

mcj 

0 

0 

1 

(2,1-12) 
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Fig,  2,1-2  Free  body  diagram  of  a  point  mass 
12 


Consider  now  an  elastomechanical  system,  such  as 
that  shown  in  Pig,  20l-3»  made  up  of  many  components 
connected  end-to-end,  each  of  which,  as  should  now  be 
apparent,  may  be  characterized  by  a  unique  transfer 
matrix o   Subscripting  the 


Flgo  2„1=3  An  elastomechanical  system 

transfer  matrices  to  reflect  the  subscripts  assigned 
to  the  state  vectors  which  describe  the  force  systems 
and  displacements  at  the  ends  of  each  component,  we  may 
write  the  following  matrix  equations; 

z2=ul,2eZl>  z3=u2,3°z2)°o0>zn=un-l,n°zn~l     (2,1-13) 

Note  that  successive  state  vectors  may  be  eliminated, 
by  substituting  for  each  its  matrix  product  equivalent,  . 
until  finally  the  matrix  relationship 

^n=Un_jL  n«°°U2  3*Ui  2°zl  remains  o   The  matrix  product, 

U'n-l  noooUo  2°u2  1  »  of  the  component  transfer  matrices 
is  itself  an  m-by-m  matrix  which  relates  the  final  to 
the  initial  state  vector  of  the  system0   It  is,  in  fact, 
the  transfer  matrix  of  the  system0 
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2a2  Elastic  Intermediate  Conditions  and  Branched  Systems 

Piping  systems  frequently  include  branches,  and 
are  supported  intermittently  along  their  length  by  hangers. 
These  branches  and  hangers  usually  have  a  very  signifi- 
cant effect  on  the  dynamic  behavior  of  a  system .   Their 
effects  may  be  considered  within  the  framework  of  the 
transfer  matrix  method e 

a0   Hangers   It  is  clear  that  the  effect  of  a  hanger  on 
a  system  is  to  introduce  discontinuities  in  the  internal 
forces  in  the  system  at  the  point  of  attachment 0 

Consider  a  spring  hanger  having  the  torsional  spring 
constants  CTX  about  the  x  axis,  and  CTY  about  the  y  axis, 
and  the  linear  spring  constant  CLX  in  the  x  direction 
For  a  Case  2  situation,  the  equations  relating  the  state 
vector  elements  on  either  side  of  the  point  of  attachment 
are: 

^R  =  ^L 

Tr  =  TL  +  CTX*S 
~WR  =-WL 

^R  «VL 
MR  =  ML  +  CTY'1(/ 


H 


or  in  matrix  notation; 


"  s 

T 

-w 

V 

= 

"y 

H 

R 

"1 

0 

0 

0 

0 

CTX 

1 
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0 

0 

U 

0 

I 

0 

0 
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0 

0 

1 

0 

0 

0 

0 

CTY 

1 

0 

0 

-CLX 

0 

0 

s 

T 

0  ¥ 
0  M, 

1 

_ 

It  should  be  apparent  at  this  point  that  a  rigid 
intermediate  support,  i0e0,  one  having  one  or  more 
infinitely  stiff  springs,  may  not  be  considered  using 
this  method .   Rigid  supports  may  be  approximated  however 
by  using  very  large  values  for  appropriate  spring  con- 
stants o   A  value  of  100,000  lbs0/ino  or  in0~lb0/  rado 
as  applicable  would  be  an  adequate  approximation  in  most 
instances,  and  it  does  not  introduce  computational  diffi- 
culties . 

b0   Branches     The  joining  of  a  branch  to  a  system 
has  an  effect  similar  to  that  of  a  spring  hanger  on  the 
force  system  in  the  main  member  at  the  point  of  joining 0 
In  fact,  it  "looks"  exactly  like  a  spring  hanger  to  the 
main  member  except,  as  will  shortly  be  seen,  its  parameters 
depend  on  frequency •   Discontinuities  are  introduced  in 
the  forces  whereas  the  main  member  deflections  on  each 
side  of  the  joint  are  the  same,  and  identical  to  those 
of  the  branch  at  the  same  point 0 


Consider  a  Cise  1  situation  and  system  shown  in 
Figo  2o2~l0 


Wain  system 


Branch 


Figo  2o2=4  Main  system  with  a  "branch 
In  the  free  body  diagram,  Figo  202-2a,  of  the  joint  of 
the  system  of  Fig0  2C2-1,  the  force  systems  acting  on 
the  joint  are  expressed  in  terms  of  their  own  coordinate 
systems „    By  a  suitable  coordinate  transformation  it 
is  possible  to  express  the  force  system  of  the  end  of 
the  branch  in  terms  of  the  main 


A       V 


MT 


Figo  2o2-2a 
member  coordinate  system . 
Fig,  202-2b. 


%\ 


^ 


\ 


Fig0  202-2b 
Such  a  situation  is  shown  in 


llhe  symbol  "  refers  to  the  fact  that  the  vector 
elements  are  expressed  in  terms  of  the  coordinates  of 
the  branch o 
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Equilibrium  considerations  and  the  orthogonality  of 
the  force  and  shear  vectors  permit  us  to  write  the 
following  equations: 

M   =  M  +  M 

R    L    B 

V  =  V  +  V 
R    L    B 

N   =  N   +  N 
R    L    B 


or  in  matrix  form 


V  + 


M 
V  = 


(2o2=0) 


Bearing  in  mind  the  continuity  of  displacements 

at  a  Joint,  it  is  clear  that  the  main  member  state 

vector  ZR  can  be  obtained  by  simply  adding  the  force 

components  of  the  branch  state  vector  Zs,  expressed 

in  the  coordinate  system  of  the  main  member,  to  the 

force  components  of  the  main  member  state  vector  Z  • 

L 

What  follows  is  the  derivation  of  a  suitable  expression 

for  the  force  components  of  the  branch  state  vector 0 

The  main  member  state  vector  Z   is  obtained  easily 

from  the  state  vector  Z  and  the  transfer  matrices  of 

the  Intervening  main  member  components 0 

At  point  C  (Figo  2  0  2-1)  there  are  three  zero  and 

1 


three  non  zero  elements  in  the  state  vector  Z 


c,„ 


■•■See  Appendix  B,  section  Bo 92  for  treatment  of  end 
conditions  which  do  not  meet  the  requirement  for  three 
zero  and  three  non  zero  elements „ 
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The  column  matrix  of  the  non  zero  components  is 
symbolized  v«   After  relating  the  branch  state  vector 
at  3  to  that  at  C  by  the  transfer  matrix  of  the  branch 
(obtained  in  the  usual  manner  from  the  transfer 
matrices  of  its  components),  and  eliminating  those 
columns  of  the  branch  transfer  matrix  which  would 
otherwise  be  multiplied  by  the  zeros  of  Zc,  we  obtain 
the  matrix  equation 

(2o2-l) 


wherein  the  state  vector  Zfi  has  been  reordered  and  par- 
titioned such  that  dLdenotes  a  sub  vector  containing 
the  elements  of  Z3  pertaining  to  deflection,  and  La 
sub  vector  containing  those  pertaining  to  force „   Expan- 
sion of  equation  202-l  gives 


aB 


and 


(2,2-2) 


(2.2-3) 


By  substituting  RJ1  °  <JR  =  v  in  equation  202»3  we 
obtain 


*B   -  R2'  Y1  °  K 


(2,2-*0 


relating  the  displacements  at  B  to  the  forces  In  the 

-1 
branch  at  B0   The  matrix  product  H^l   is  calleci  a 

spring  matrix  because  of  the  function  that  it  serves 0 
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although  any  consistant  system  will  do,  it  is  con- 
venient here  to  establish  a  particular  method  for  measuring 
the  angle  at  which  a  branch  joins  the  main  member «   The 
one  defined  is  consistant  with  that  used  in  VIPIPEo 
Accordingly,  -the  angle  at  which  a  branch  joins  the  main 
member  is  measured  counter-clockwise  from  the  x  axis  of 
the  main  member  as  it  is  oriented  immediately  before 
the  intersection,  to  the  x  axis  of  the  branch  as  it  is 
oriented  at  the  point  of  intersection.   The  angle  0  in 
Figo  2o2-l  is  drawn  in  conformance  with  this  definition0 

From  the  matrix  relationships 

dB  =  G^  dB 
and 

pb  =  v  h 

where  for  a  planar  piping  system 


Gi  = 


cos0   -sin0   0 

sin0   cos0   0 

0      0    1 


and  G0   =■ 


1   0 
0  cos0 
0   sin0 


0 

-sin0 
cos0 


the  spring  matrix  (symbolized  S) 
S  =  G2°  R2°  R1"  .  G1 

is  obtained  which  relates  the  force  and  displacement 
elements  of  the  branch  transfer  matrix  Zfi  (expressed  in 
terms  of  the  coordinates  of  the  main  member) 0   Thus 
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It  is  now  a  simple  matter  to  construct  the  transfer 
matrix  relating  the  main  member  state  vectors  on  either 
side  of  the  joint •   Recall  from  equation  2o2-0  that 

?R  =  ?L  +  PB  =  PL  +  3o*3 
since 

dD  =  dT  =  dp 


AR 


"B  , 


i ;  o 

S  !  I 


ur  zL 


where  we  use  the  symbol  U  to  denote  the  point  matrix 
which  has  just  been  derived . 

Thus  a  branch  is  treated  as  a  point  transfer  matrix 
in  constructing  the  transfer  matrix  of  the  systemo 

2»3  Mode  Frequencies  and  the  Frequency  Determinant 

A  principal  application  of  the  method  of  transfer 
matrices,  and  the  one  with  which  we  are  concerned,  is 
to  the  vibration  analysis  of  elastomechanical  systems 
for  the  purpose  of  obtaining  their  natural  frequencies » 
When  transfer  matrices  are  constructed  for  systems 
having  mass,  certain  of  their  elements  are  frequency 
dependent  (recall  the  transfer  matrix  of  a  point  mass 
derived  earlier) «   Proper  application  of  the  system 
boundary  conditions  to  these  transfer  matrices  yields 
a  matrix  of  order  n'2  (from  a  transfer  matrix  of  order  n) 
which  can  be  used  to  obtain  the  natural  frequencies  of  a 
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system 0  Its  determinant  Is  termed  the  frequency  determi- 
nant of  the  system.  The  concept  is  best  described  by 
an  example • 

Consider  a  case  2  situation  wherein  a  system  has 
a  system  transfer  matrix  \J3   and  boundary  conditions 
embodied  in  its  left  and  right  end  state  vectors  ZL  and 
ZRo   Assume  that  the  left  end  is  fixed  and  that  the  right 
end  is  free»   The  matrix  equation  relating  the  left  and 
right  end  state  vectors  is 


5 

0 

-w 

o 
ol 


ul,l 


ul,6 


'JH  Lu6,l    oooo  ii6,6 

from  which  may  be   constructed   the  matrix  relationship 

0 


u2,2  "2,5  u2,6 
u5,2  "5,5  u5,6 
u6,2     »6,5     u6,6 


(2o3-D 


Clearly,  a  non-trivial  solution  requires  that  the  third 
order  submatrix  of  equation  203~1  must  vanish 0   The 
determinant  of  the  submatrix  is  the  system  frequency 
determinant  mentioned  above „   The  frequencies  for  which 
its  value  is  zero  are  the  natural  frequencies  of  the 
system  described  by  the  transfer  matrix  Ug  and  the 
boundary  condition  state  vectors,  ZL  and  ZR,  from  which 
2  1 


it  was  derived  a 


2 ,k     The  Digital  Computer  and  the  Method  of  Transfer  Matrices 

Even   for  a  relatively  simple  system,  the  elements 
of  the  frequency  determinant  become  unwieldly  polynomials 
in  radian  frequency 0   Determination  of  the  roots  of  the 
equation  resulting  from  an  expansion  of  the  frequency 
determinant  is  usually  quite  impractical  and  recourse 
must  be  made  to  a  numerical  solution o   In  a  numerical 
solution  the  frequency  determinant  is  constructed  and 
evaluated  for  selected  frequencies  over  a  range  of 
interest.   Its  value  is  plotted  as  a  function  of  fre- 
quency and  the  natural  frequencies  of  the  system  which 
it  represents  are  the  zeros  of  the  plot  (Fig0  203-U° 


A    (sihc  value  of 


Flg-„2o3~l  Frequency  determinant  as  a  function  of  frequency 

Use  of  a  digital  computer  permits  rapid  and  accu- 
rate analysis  of  complex  systems 0   Program  VIPIPE 
presented  in  Appendices  A  through  C,  for  instance,  will 
accept  for  analysis  planar  piping  systems  having  up  to 
fifty  components  including  as  many  as  twenty  branches 
and  spring  hangers „   It  treats  distributed  mass  as  such 
and  considers  the  effects  of  shear  deflection  and  rota- 
tional inertia .   In  a  typical  analysis,  the  first  five 
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natural  frequencies  were  found  for  a  system  having 
thirteen  components,  including  two  spring  hangers  and 
two  branches,  at  the  cost  of  five  minutes  of  machine 
time  on  a  CDC  160*1-  computer . 

It  is  through  the  use  of  high  speed  digital 
computers,  that  the  full  potential  of  the  method  of 
transfer  matrices  in  this  application  as  well  as  in 
otherslmay  be  realized e 


••■See  ref  o  1  for  a  comprehensive  treatment  of 
applications  of  the  method  of  transfer  matrices  to 
other  problems  in  elastomechanics * 
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CHAPTER  III 


DISCUSSION 


3.1   Accuracy  of  Method 

Confidence  in  the  accuracy  of  the  method  was 
established  by  using  it  to  determine  the  natural 
frequencies  of  several  systems  for  which  these  fre- 
quencies could  be  obtained  by  other  means,  i0eo,  by 
means  of  either  a  closed  form  analytical  solution  or 
by  recourse  to  the  literature  [2]  o   The  average  dif- 
ference between  system  natural  frequencies  obtained 
using  VTPIPE  and  those  obtained  by  analytical  means 
was  found  to  be  less  than  O0064$,  while  the  average 
difference  between  system  natural  frequencies  obtained 
using  VIPIPE  and  those  obtained  from  the  literature 
was  found  to  be  less  than  05$» 

Details  of  the  systems  analyzed  and  the  results 
obtained  may  be  found  in  Appendix  E,  section  Eo210   The 
analyzed  systems  consisted  of  single  spans  of  straight 
pipe  having  various  combinations  of  boundary  conditions . 

The  published  papers  of  E0  Volt era  [V]  give 
methods  for  determining  the  natural  frequencies  of 
curved  filaments,  but  the  methods  are  sufficiently 
involved  as  to  make  their  application  in  obtaining 
some  comparison  frequencies  a  major  undertakingo 
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Time  precluded  assumption  of  the  task.   The  literature 
[2,3]  does  provide  means  for  computing  the  fundamental 
frequency  of  partial  rings  with  fixed  ends0  The  dif- 
ference between  the  fundamental  frequencies  determined 
using  VIPIPE  and  those  computed  using  the  method  of 
the  literature,  was  found  to  be  less  than. 3$  for  in~ 
plane,  and  2„3$  for  out-of-plane  vibrations,1  Details 
of  the  system  analyzed  and  the  results  obtained  may  be 
found  in  Appendix  E,  section  E02^0 
3.2   Integrity  of  Solution 

Confidence  in  the  integrity  of  the  method  was 
established  by  conduct ing"mirror  image"  analyses  of 
almost  every  configuration  considered „   Construction 
of  the  transfer  matrix  of  a  system  is  directional  in 
nature,  in  that  an  end  of  the  system  is  chosen  as  a 
starting  point,  and  the  transfer  matrix  of  the  system 
is  developed  by  premultiplying  the  incomplete  system 
transfer  matrix  by  the  transfer  matrix  of  each  system 
component  in  the  order  that  each  component  is  encoun- 
tered in  tracing  the  system  from  that  starting  point • 
It  follows  therefore  that  an  error  in  the  structure  of 
any  of  the  component  transfer  matrices,  or  in  the  method 
or  order  of  their  combination  in  constructing  the 
system  transfer  matrix,  would  reveal  itself  in  the  form 


See  footnote  to  section  E«2ty  of  Appendix  E. 
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of  inconsistent  results  if  the  system  were  analyzed 
twice,  beginning  the  first  time  from  one  end,  the  second 
time  from  the  other  end0   Such  a  dual  analysis  is 
referred  to  herein  as  a  "mirror  image"  checko 

In  general,  the  disagreement  between  natural 
frequencies  obtained  for  a  system  and  for  its  mirror 
image  xvas  found  to  be  less  than  o0003/£» 

An  added  check  on  solution  integrity  was  obtained 
by  analyzing  a  single  branch  system  treating  first 
one  member  then  another  as  the  branch 0   The  average 
difference  here  was  on  the  order  of  O027$  for  in-plane 
vibration,  and  o00?^  for  out-of-plane  vibrationc 

Details  of  the  systems  analyzed  and  of  the  results 
obtained  may  be  found  in  Appendix  E,  sections  E031  and  E„32. 
3.3  Computer  Imposed  Limitations 

The  number  of  natural  frequencies  that  may  be  found 
for  a  given  system  using  VTPIPE,  is  a  function  of  the 
floating  point  significant  figure  capacity  of  the  com- 
puter usedo   The  significant  figure  capacity  required 
to  solve  for  a  given  number  of  natural  frequencies  of 
any  system  is  a  function  of  the  characteristics  the 
system;  and  for  all  systems,  the  significant  figure 
capacity  required  increases  with  the  number  of  the  mode 
sought ■   Thus,  for  example,  for  a  given  system,  the 
number  of  significant  figures  required  in  order  to 
determine  the  third  mode  frequency  is  greater  than 
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that  required  to  determine  the  first  .node  frequency, 
and  while  the  first  ten  natural  frequencies  of  one 
system  may  be  obtained  without  difficulty,  the  signifi- 
cant figure  capacity  of  the  computer  may  be  inadequate 
to  determine  the  fourth  mode  frequency  of  another0 
Program  VIPIPE  is  designed  to  terminate  a  solution 
when  the  number  of  significant  figures  required  for  its 
continuance  exceeds  that  available  in  the  machine ■ 

The  use  of  multiple  precision  arithmetic  would 
extend  VIPIPE »S  utility  considerably 0 
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APPENDIX  A 


DESCRIPTION  UF  PiIOGRaM 


Ad  General  Remarks 

Program  VIPIPE  is  a  CDC  FORTRAN  60  [$]  language 
digital  computer  program  designed  for  use  in  the  vibra- 
tion analysis  of  planar  piping  systems 0   In  and  out-of- 
plane  undamped  natural  frequencies  may  be  found  through 
its  use»   The  analysis  is  performed  by  an  application 
of  the  method  of  transfer  matrices  wherein  each  component 
of  a  system  is  characterized  by  a  matrix  called  its 
transfer  matrix •   The  entire  system  is  characterized 
by  the  system  transfer  matrix  obtained  through  an 
ordered  multiplication  of  the  component  transfer 
matrices.   Considerable  flexibility  with  regard  to 
the  mathematical  model  employed  to  describe  the  system, 
and  the  nature  of  simplifying  assumptions  accepted  in 
its  vibration  analysis,  is  provided  in  the  program0 
For  instance,  although  the  method  of  solution  yielding 
the  greatest  accuracy  coupled  with  least  machine  time 
combines  a  distributed  mass  treatment  of  straight 
sections  with  a  lumped  mass  treatment  of  curved  sections, 
and  considers  the  effects  of  both  shear  deflection  and 
rotational  inertia,  the  user  may  elect  to  use  a  lumped 
mass  treatment  of  all  sections  and  to  consider  or  neglect, 
singly  or  together,  the  effects  of  shear  deflection  and 

rotational  inertia 0 
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Problem  solution  proceeds  by  an  iterative  process 
with  radian  frequency  the  independent  variable*   The 
sign  of  the  frequency  determinant  •'-is  used  to  control 
the  search  for  and  convergence  to  the  system  natural 
frequencies,,   The  frequency  increment  is  computed 
internally,  as  is  the  solution  acceptability  criterion 
Both  are  accessible,  however,  and  may  be  controlled  by 
the  user*   Subsectioning  of  components  treated  by  lumped 
mass  methods  is  carried  out  internally  as  a  function  of 
component  L/D  (length-to-diameter)  ratio  and  the  highest 
mode  sought o 

The  program  has  the  capability  to  analyze  piping 
systems  comprised  of  a  maximum  of  fifty  components,  and 
having  as  many  as  twenty  branches  and  twenty  hangers  0 

A  component,  for  program  purposes,  is  defined 
as  a  curved  section,  i0eo,  a  section  of  constant  curva- 
ture, a  straight  section,  a  spring  hanger,  a  flange  or 
a  coupling 8 ^  Although  straight  and  curved  sections, 
flanges  and  couplings  must  each  be  of  constant  outside 
diameter  and  wall  thickness,  these  dimensions  may  vary 
from  component  to  component s   Similarly,  intensive 
properties,  l,e«,  shear  and  elastic  moduli,  and  density, 
may  not  vary  within  a  component,  but  may  vary  from  com- 
ponent to  component . 


1 

See  Chapc  2,  2,3  of  the  main  body0 
^See  App0  B,  section  B09^0 
^See  Appc  B,  section  B<>95o 
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A  shear  area  form  factor  of  one  half  is  built  into 
the  program*,   It  may  be  altered  at  the  user's  discretion 
by  altering  one  carde1 

Five  common  end  or  boundary  conditions,  for  the 
ends  of  the  main  member  and  the  remote  ends  of  any 
branches,  may  be  specified  by  a  simple  coded  input e 

Unusual  boundary  conditions  may  be  considered  by  recourse 

2 

to  a  more  complicated  input 0 

Systems  having  rigid  intermediate  supports  and/or 
non  linear  elements  may  not  be  analyzed  by  the  use  of 
this  program,,   The  number  of  natural  frequencies  found 
for  a  given  system  is  a  function  of  the  significant 
figure  capacity  of  the  computer  employed  since  the  num- 
ber of  significant  figures  required  in  a  solution  is  a 
function  of  the  characteristics  of  the  system  under 
analysis,  and  for  any  system  increases  with  increasing 
mode  number o 
Ao2  Program  Structure 

The  program  consists  of  a  six  section  main  body  and 
twenty-four  subroutines „   If  sine*  cosine,  sinh,  and 
cosh  function  routines  are  not  available  in  the  FORTRAN 
library  they  must  be  provided  by  other  means  since, 
although  required  by  the  program,  they  are  not  explicitly 
contained  therein •   The  function  of  each  section  of  the 
main  body  and  of  each  subroutine  is  given  below8 


1See  Appc  B,  section  B091« 
2See  App.  B,  section  B„64o 
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A  ,21  Main  Body  Sections  and  their  Functions 

INPUT       Control  read-in  of  data,  allocate  storage 
locations  for  arrays  and  subscripted  vari- 
ables 0 

INVARIANTS  Compute,  and  assign  appropriate  subscripts 
to,  such  system  component  characteristics 
as  moment  of  inertia,  and  radius  of  gyra- 
tion, 

CONTROL  1   Construct  transfer  matrix  of  the  system 
for  the  in-plane^case  by  calling  for 
various  subroutines  in  a  sequence  governed 
by  the  system  geometry,  the  mathematical 
model  specified,  and  the  desired  assumptions,, 

CONTROL  2   CONTROL  2  performs  the  same  function  as 
CONTROL  1  but  for  the  out-of~plane2casec 

ITERATION   Apply  boundary  conditions  to  the  system 

transfer  matrix  constructed  by  CONTROL  1  or 2. 
so  as  to  obtain  the  system  frequency 
determinant o   Evaluate  the  determinant 
and  utilize  its  sign-^in  controlling  the 
iterative  search  for  mode  frequencies. 
Initiate  output  when  the  specified  fre- 
quency range  has  been  fully  explored,  the 
required  number  of  modes  found,  or  the 
significant  figure  limit  of  the  com outer 
reached » 

OUTPUT      Control  output  format «, 

A o 2 2   Subroutines  and  their  Fun c t Ions 


SUB3EC 


DI6TM 


Subdivide  components  to  be  treated  by- 
lumping  mass  on  the  basis  of  L/D   ratio 
(or  ratio  of  radius  of  curvature  to  diam- 
eter ^and  included  angle  of  arc,  as 
applicable. ) 

Construct  transfer  matrix  of  straight' 
sections  having  distributed  mass  for  Case  10 


Henceforth  Case  1 „ 
-Henceforth  Case  2o 
3see  App0  B,  section  B09^o 
^Henceforth  RHO/D  where  RHO  is  the  radius  of  curvature 

of  the  component u 
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Dlori-IO      Construct  transfer  matrix  of  straight 
sections  having  distributed  mass  for 
Case  2o 

SFIELD      Construct  the  transfer  matrix  of  straight 
nassless  subsections  for  Case  1„ 

3FIELG      Construct  the  transfer  matrix  of  straight 
massless  subsections  for  Case  2, 

CFIELD      Construct  the  transfer  matrix  of  curved 
massless  subsections  for  Case  1. 

CFIELO      Construct  the  transfer  matrix  of  massless 
curved  subsections  for  Case  20 

POINT  Construct  the  transfer  matrix  of  a  point 
mass  for  Case  10 

PuINO  Construct  the  transfer  matrix  of  a  point 
mass  for  Case  2* 

MATMUL      Construct  the  transfer  matrix  of  those 

components  treated  by  lumping  mass  through 
the  ordered  matrix  multiplication  of  the 
transfer  matrices  of  massless  subsections 
and  point  masses 0 

RIGID  Construct  the  transfer  matrix  of  a  rigid 
straight  section  having  distributed  mass 
for  Case  1« 

RIGIO  Construct  the  transfer  matrix  of  a  rigid 
straight  section  having  distributed  mass 
for  Case  20 

3TIFC0      Construct  the  transfer  matrix  of  a  massless 
rigid  corner  (essentially  a  transformation 
of  coordinates  at  an  elbow)  for  Case  10 

STIFOO      Construct  the  transfer  matrix  of  a  massless 
rigid  corner  (essentially  a  transformation 
of  coordinates  at  an  elbow)  for  Case  2« 

HANGER      Construct  the  transfer  matrix  of  a  spring 
hanger  for  Case  1 « 

HANGEO      Construct  the  transfer  matrix  of  a  spring 
hanger  for  Case  20 

BRANCH      Construct  a  transfer  matrix  from  the  trans- 
fer matrix  of  a  branch  and  its  remote  end 
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boundary  conditions,  which  permits  a 
branch  to  be  treated  as  a  component  of 
the  main  member y  for  Case  lc 

BRaKCO      Construct  a  transfer  matrix  from  the 
transfer  matrix  of  a  branch  and  its 
remote  end  boundary  conditions,  which 
permits  a  branch  to  be  treated  as  a  com- 
ponent of  the  main  member  for  Case  2<, 

STAVEC      Construct  boundary  condition  vectors 
from  coded  input  for  Case  10 

STAVEO      Construct  boundary  condition  vectors 
from  coded  input  for  Case  20 

INVERT      Invert  a  matrix 0 

FINMAT      Construct  the  transfer  matrix  of  the 

system  by  successive  premultiplication 
of  the  partial  system  transfer  matrix  by 
the  transfer  matrix  of  each  system  com- 
ponent as  soon  as  each  of  the  latter  is 
constructed « 

FINBRA      Construct  the  transfer  matrix  of  a  branch 
by  successive  premultiplication  of  the 
partial  branch  transfer  matrix  by  the 
transfer  matrix  of  each  branch  component 
as  soon  as  each  of  the  latter  is  constructed., 

DELMA       Construct  the  frequency  determinant  of 

the  system  by  applying  the  boundary  con- 
ditions of  the  main  member  to  the  com- 
pleted system  transfer  matrix 0 


Ao3   Program  Nomenclature 
A. 31   Main  Body 
A 


AA 


AAriU 


Array /matrix o   Transfer  matrix  of  a 
massless  subsection0 

Discriminant o   A  control  in  the  conver- 
gence sequence  <• 

In  INPUT;  density  (lbs./cUoft „ )  of  the 
first  system  component 0 
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In  INVARIANTS;  average  mass  per  unit 
length  of  the  main  member! 

AMU         In  INPUT |  component  density  (lbs/cu<,  f 1 0 )  . 
In  remainder  of  program;  component  mass 
per  unit  lengths 

iiD         Average  outside  diameter  of  the  main 
member u2 

iUJUii        Current  computed  iteration  frequency- 
increment  o 

APR        Average  moment  of  inertia  of  the  main 
member o 3 

AE         Elastic  modulus  of  the  first  system 
component o 

nG  Shear  modulus  of  the  first  system  com- 

ponent e 

AID         Average  inside  diameter  of  the  main  member^ 

AIX         Component  radius  of  gyration  about  the 
x  axiSo 

kll  Component  radius  of  gyration  about  the 

y  or  z  axiSo 

AJ  Component  moment  of  inertia  about  the 

y  or  z  axiSo 

AJT         Component  polar  moment  of  inertia e 

AL         Total  length  of  main  member 0 

AMYK        Discriminant o   Controls  choice  of  mathe- 
matical model,  i0ec,  lumped  or  distributed 
masso 

B  MatriXo   Transfer  matrix  of  a  point  mass0 

BB         Discriminant „   Controls  Iteration  fre- 
quency during  convergence 0 


solution  acceptability  criteriom 
2~^Employed  in  computing  the  initial  frequency  increment 
and  solution  acceptability  criterion0 
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3D0M  Initial  frequency  increment o 

BAMU  Component  density  (lba/cu,  ft0) 

BM  Number  of  last  mode  found c 

dR  Initial  solution  acceptability  criterion, 

CC         Discriminant ,   Designates  final  iteration 
during  convergence 0 

CLA  Linear  spring  constant  of  a  hanger  in 

the  x  direction.,1 

CLY         Linear  spring  constant  of  a  hanger  in 
the  y  direction, 

CLZ         Linear  spring  constant  of  a  hanger  in 
the  z  dlrection,3 

CR  Current  value  of  the  solution  acceptability 

criterion, 

CTX         Torsional  spring  constant  of  a  hanger  about 
the  x  axis0* 

CTY         Torsional  spring  constant  of  a  hanger  about 
the  y  axise5 

CTZ         Torsional  spring  constant  of  a  hanger  about 
the  z  axis,° 

D  Component  outside  diameter, 

DD  Discriminant,   Controls  read-in  and  stor- 

age of  the  elastic  and  shear  moduli  and 
density  of  components • 

DI  In  INPUT;  component  wall  thickness, 

in  remainder  of  program;  component  inside 
diameter, 

DT  Component  wall  thickness, 

DOM         Current  frequency  increment  during  con- 
vergence. 


1~"No  coupling  is  permitted  in  action  of  hanger.   For 
treatment  of  hangers  in  which  coupling  occurs  see 
App.  B,  3,9. 


£  Component  elastic  modulus , 

£E         Average  elastic  modulus  of  the  main 
member. 

FF.dC        Shear  area  form  factor  0 

FR         First  mode  frequency  of  a  fictitious 

straight  beam  having  fixed  ends,  physical 
properties  and  transverse  dimensions 
equal  to  an  average  of  the  characteristics 
of  the  main  member  to  which  it  is  equal 
in  length*, 

G  Component  shear  modulus 0 

GG         Discriminant ,   Control  in  test  for  approach 
to  significant  figure  limit  of  computer,, 

HDEL        Reference  value  of  the  frequency  deter- 
minant (significant  figure  limit  test). 

HM         Number  of  highest  mode  sought, 

I  Subscript c 

II  Subscript  which  associates  a  branch  with 
the  angle  at  which  it  joins  the  main 
member. 

IOP         Discriminant o   Controls  analysis  to  pro- 
vide Case  1  and/or  Case  2  mode  frequencies 
as  specified. 

JK         Discriminant.   Controls  solution  sequence 
during  Case  1  or  2  analysis. 

KK         Discriminant.   Controls  output  to  provide 
a  print-out  that  an  analysis  has  been 
terminated  due  to  significant  figure  limi- 
tations when  appropriate. 

L  Iteration  number. 

LM  Number  of  current  iteration. 

M  Number  of  mode  currently  being  sought. 

MH  Number  of  highest  mode  found. 

N  Subscript  which  associates  a  branch  with 

its  remote  end  boundary  condition. 
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NBR         Number  of  branches  in  a  system. 

M         Discriminant o   Controls  output „ 

NMBR        Total  number  of  system  boundary  conditions. 

1S1N3R        Discriminant*,   Assigns  a  component  to 
the  main  member  or  to  a  branch,  and  if 
the  latter,  indicates  its  position  in 
the  branch o 

IMS  Total  number  of  system  components  . 

OMEGA       Current  iteration  frequency 0 

OiiEGAC      Convergence  reference  frequency. 

OiflEGAG      Upper  limit  of  frequency  range  of  investi- 
gation » 

QMEGAI      Investigation  starting  frequency. 

OMBGAM      Mode  frequency. 

P  Sum  of  the  numerical  values  of  the  spring 

constants  of  a  component. 

PDEL        Value  of  the  frequency  determinant  used 
as  a  reference  during  iteration. 

PHI         Angle  at  which  a  branch  joins  the  main 
member. 

R  Reciprocal  of  component  flexural  rigidity. 

RRO         Component  radius  of  curvature. 

RI  Discriminant.   Controls  consideration 

of  rotational  inertia. 

3D  Discriminant.   Controls  consideration  of 

shear  deflection. 

3V  Coded  boundary  condition  for  Case  1. 

SVTN  Non-coded  boundary  condition  for  Case  1. 

SVO  Coded  boundary  condition  for  Case  2. 

SVOP  Non-coded  boundary  condition  for  Case  2. 

TDEL        Discriminant.   Controls  significant  figure 
cut  off. 
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THETA 

THETR 

TL 
U 

UU 
V 

W 
Z 

A.  32  Subroutines 
a.  SUBSEC 
RHOD 

RA 
RD 


Included  angle  of  arc  of  curved  com- 
ponents in  degrees  in  input,  radian 
measure  in  remainder  of  program. 

Included  angle  of  arc  of  curved  com- 
ponents in  degrees o 

Total  length  of  a  component 0 

Array  name.   Component  transfer  matrix • 

Array  naie.   System  transfer  matrix 0 

Array  name.   System  or  branch  transfer 
matrix. 

Array  name0   Branch  transfer  matrix,, 

Number  of  subsections  of  a  component. 


Absolute  value  of  radius  of  curvature  of 
component o 


Component  L/D  ratio • 
Component  RHO/D  ratio e 


b.  DITM  and  DISTMO 


ARA 


AREA 

ALl^ 

AL2 

CBS 

CHL1 

CL2 

COO 

C01 

C02 


Component  cross  sectional  area. 
Component  equivalent  shear  area. 


C03 

BE 

CL1 

s 

SBE 

SHL1 

SL2 

SPT 

W 


> 


Array  elements  and/or  computed  constants 
used  in  evaluating  array  elements. 
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}R       Component  flexural  rigidity. 

W         Component  torsional  stiffness  (DISTMO). 

3 FIELD  and  5FIEL0 


ARa       Component  cross  sectional  area  (SFIELD  only). 
AREA      Component  equivalent  shear  area. 
DL        Length  of  subsection,, 

CFIELD  and  CFIELO 


RHOD 

SD 

V 


Component  cross  sectional  area. 

Component  equivalent  shear  area  (CFIELD  only). 


Array  elements  and/or  constants  used 
in  computing  array  elements e 


Component  shear  stiffness  (CFIELO  only). 

Included  angle  of  arc  of  subsection0 

Absolute  value  of  radius  of  curvature 
(CFIELD  only)o 

Component  longitudinal  stiffness  (CFIELD 
only) . 

Signed  component  radius  of  curvature 
(CFIELO  only). 

Discriminant .   Controls  consideration 
of  shear  deflection  (CFIELO  only)0 

Factor,   Imparts  sign  to  certain  array 
elements. 


e.  POINT,  POINQ.  PIGID.  and  PIGIO 

AM        Total  mass  of  subsection  or  component. 

f.  STIFCQ  and  STIFOO 


Factor,,   Imparts  sign  to  certain  array 
elements. 
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g.   BRANCH  and  BRANCO 

W        Branch  transfer  matrix  (3HAKCH) 

UV        Branch  transfer  matrix  (BRAN CO), 

3VB       Array  name.   Remote  end  boundary  condi- 
tion vector  (BRANCH), 

3R       Array  name ,   Remote  end  boundary  condi- 
tion vector  (BRANCO). 

Rl ,  R2    Array  names,   Submatrices  of  branch 
transfer  matrix, 

R  Array  name.  Matrix  product  of  R2  °  Rl"1 
in  first  usage,  spring  matrix  of  branch  in 
second  usage  in  BRANCH • 

S         Array  name.   Spring  matrix  of  branch 

BRANCO, 

Gl ,  G2    Coordinate  transformation  matrix  elements, 

U  Array  name.  Transfer  matrix  which  repre- 
sents a  branch  as  a  component  of  the  main 
member. 
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APPENDIX  B 
INSTRUCTIONS  FOR  PROGRAM  USE 

B.l   General  Remarks 

Appendix  B  is  intended  to  provide  complete  instruc- 
tions for  using  program  VIPIPEo  To  that  end,  since  for 
problem  programming  purposes,  most  system  components  must 
be  identified  with  respect  to  a  coordinate  reference,  a 
description  of  the  reference  system  employed  is  given  first. 
It  is  followed  by  a  summary  of  data  required,  a  description 
of  the  component  identification  system,  an  explanation  of 
the  nature  of  a  state  vector  in  terms  of  whose  elements 
the  system's  boundary  conditions,  and  its  displacements 
and  internal  force  systems  are  described,  and  a  description 
of  the  data  input  format.   The  appendix  concludes  with 
some  comments  on  the  options  available  in  the  program. 

For  the  purpose  of  setting  up  a  problem  for  solution, 
the  longest  continuous  run  of  pipe,  or  the  one  from  which 
all  branches  emanate  if  these  are  not  one  and  the  same, 
must  be  designated  as  the  main  member.   A  system  component 
is  defined  as  a  section  of  pipe,  straight  or  of  constant 
curvature,  an  elbow,  a  spring  hanger,  a  coupling-^or  a 
flange.   Diameter,  wall  thickness,  elastic  and  shear  moduli, 
and  density  may  vary  from  component  to  component,  but  may 
not  vary  within  a  component.   Since  both  torsional  and 
linear  spring  rates  (for  hangers)  appear  internally  as 


^See  App.  B,  section  B.95. 
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elements  of  a  matrix  used  in  obtaining  matrix  products, 
they  must  have  finite  values 0   Thus  perfectly  rigid 
hangers  may  not  be  treated  exactly,  but  they  can  be 
approximated  by  using  very  large  spring  constants . 
B.2   Coordinate  System 

In  accordance  with  the  procedures  given  below,  a 
local  coordinate  system  is  constructed  at  each  point  of 
interest  in  the  system,,   These  points  are  usually  the 
ends  of  each  component,  the  beginning  and  end  of  the  main 
member,  and  both  ends  of  each  branch »   The  coordinate 
systems  are  established  with  the  piping  system  quiescent. 
Once  established,  they  remain  fixed  in  space . 

In  establishing  a  coordinate  system,  one  end  of  the 
main  member,  henceforth  called  the  left  end,  must  be  chosen 
as  a  starting  point «. 

A  right-handed  coordinate  system  is  used  whose  x 
axis,  always  directed  away  from  the  left  end,  is  the  local 
tangent  to  the  central  axis  of  the  main  member «   The  origin 
of  the  coordinate  system  lies  on  the  main  member  central 
axis,  and  the  xy  plane  coincides  with  the  quiescent  plane 
of  the  system,,   The  z  axis  is  always  directed  in  the  same 
direction,  and  this  direction  shall  be  called  downward. 
The  direction  of  the  local  y  axis  is  chosen  so  as  to  com- 
plete the  right-handed  system,, 

central  axis 


Fig0  Bo2   Coordinate  System 
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The  starting;  point,  for  reference  purposes,  of  a 
branch  is  taken  as  its  far  end,  i0e0,  the  end  farthest 
removed  from  the  intersection  of  the  branch  vjith  the  main 
member.   The  plane  of  the  branch  coincides  with  the  xy  plane „ 
The  positive  z  axis  of  the  branch  and  main  member  are 
parallel  and  have  the  same  sense B   The  y  axis  is  oriented 
so  as  to  complete  a  right-handed  system0 
B.3   Branch  Joining  Angle 

The  angle  at  which  a  branch  joins  the  main  member  is 

measured  counter-clockwise  in  degrees  from  the  x  axis  of 

the  main  member  as  it  is  oriented  immediately  before  the 

intersection,  to  the  x  axis  of  the  branch  as  it  is 

oriented  at  the  point  of  intersection0 

B.4  Information  Required  to  Program  Problem 

B.*fl  General 

a0   number  of  highest  mode  sought. 

b0   frequency  range  over  which  investigation  is 

to  be  performed  (initial  and  final  frequencies 
in  radians  per  second) • 

c.  whether  in-plane,  out-of -plane  or  both  in- and 
out-of-plane  mode  frequencies  are  sought 0 

d.  whether  employment  of  a  distributed  or  of  a 
lumped  mass  model  is  desired „ 

e.  whether  it  is  desired  to  consider  or  ignore 
singly  or  together  the  effects  of  shear 
deflection  and  rotational  inertia . 

f o   boundary  or  support  conditions  at  the  ends 
of  the  main  member  and  at  the  remote  end  of 
each  branch . 

go   number  of  components  in  the  system0 

ho   number  of  branches  in  the  system0 
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B.^2   All  Components  axcent  Hangers 

a„   outsi.de  diameter  ( Inches )0 

bo   wall  thickness  ( inches )0 

Co   density  (lbsu/cu0  f t>)  <> 

do   shear  modulus  (psi)0 

e.   elastic  modulus  (psi)0 
Bo^-3  Straight  Sections 

a0   length  (inches) 0 

B.^4  Curbed  Sections 

a0   signed  radius  of  curvature  (negative  if 
center  of  curvature  lies  on  the  negative 
y  axis,  positive  if  on  the  positive  y  axis ) • 

b.   arc  central  angle  (degrees )0 

B  o  4-5   Hangers1 

a0   linear  spring  constant  in  x  direction 

(lbso/injo 

bo   linear  spring  constant  in  y  direction 
(lbSo/ln0)o 

Co   linear  spring  constant  in  z  direction 

(lbs./in,), 

do   torsional  spring  constant  about  x  axis 
(in0~lbSo/rad0)o 

e,   torsional  spring  constant  about  y  axis 
(ine~lbso/radc)o 

f0   torsional  spring  constant  about  z  axis 
(in0-lbso/rado)o 

3,5   Component  Identification  System 

Components  are  identified  by  two  numbers „   The  first, 

a  component  sequence  number,  is  used  to  order  data  cards 

containing  component  extensive  and  Intensive  property  data. 


!See  App.  B„,  section  B096,  ChapQ  2,  section  2Q2a0 
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The  second,  a  branch  component  Identification  number,  is 
read  in  as  data  and  serves  to  indicate  whether  or  not  a 
component  is  a  member  of  a  branch  and  if  so,  the  position 
it  occupies  in  the  branch 0 

In  addition,  all  branches,  if  any,  are  assigned  a 
sequence  number  which  is  used  to  order  data  and/or  data 
cards  pertaining  to  branch  boundary  conditions  and  joining 
angles. 

The  methods  of  assigning  these  numbers  are  given 
below o 
B«51  Component  Sequence  Number 

All  system  components  bear  a  component  sequence 
number .   The  number  one  is  assigned  to  the  leftmost  com- 
ponent of  the  main  member ,   Each  succeeding  component 
encountered  in  tracing  the  main  member  is  assigned  a 
number  one  greater  than  that  assigned  to  its  lefthand 
neighbor •   If  and  when  a  branch  is  encountered  its  com- 
ponents are  also  assigned  sequence  numbers  beginning  with 
the  component  farthest  removed  from  intersection  of  the 
branch  with  the  main  member o   That  component  is  assigned 
a  number  one  greater  than  that  assigned  to  the  last  com- 
ponent of  the  main  member  encountered  prior  to,  and  abut- 
ting on,  the  junction.,   when  more  than  one  branch  joins 
the  main  member  at  the  same  point,  the  user  is  free  to 
choose  the  order  in  which  desired  to  take  up  the  branches 
for  the  purpose  of  assigning  sequence  numbers  to  their 
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components.   Having  established  the  order,  the  first 
branch  is  treated  as  If  it  were  the  only  branch 0   The 
remaining  branches  are  treated  in  exactly  the  same  manner 
except  that  the  component  of  the  preceeding  branch  having 
the  highest  sequence  number,  i0e0,  the  one  abutting  the 
joint,  assumes  the  role  performed  by  the  last  encountered 
main  member  component  in  treating  the  first  branch  (see 
Figo  B.5D. 

-Left  end 
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-f   TT^TO 


3-f  ©-spring  hanger 


Figo  B51   Component  Sequence  Numbers 
B„52   Branch  Component  Identification  Number 

Assignment  of  a  component  to  a  branch  is  indicated 
by  the  value  of  a  fixed  point  constant  (NNBR)  coded  as  follows: 


NN3R=1 


NNBR=2 


c .   NN3R=3 


if  component  is  first  member  of  a 
branch,  i0e0,  the  component 
farthest  removed  from  the  inter- 
section0 

if  component  is  neither  the  first 
nor  the  last  member  of  the  branch, 

if  the  component  is  the  last 
member  of  the  branch . 
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d0   NNBR=-1    If  the  branch  has  only  one  com- 
ponent , 

ee   IviWBR-0     all  main  member  components., 

Bo 53  Branch  Sequence  Number 

Branches  are  numbered  sequentially  as  they  are 

encountered  In  tracing  the  main  member 0   If  two  or  more 

branches  join  the  main  member  at  the  same  point,  the  branch 

whose  components  have  been  assigned  the  lowest  sequence 

numbers  is  assigned  the  lower  or  lowest  branch  sequence 

number  (see  Figo  3c53)° 
_0 0 


Fig*  B053  Branch  component  identification 
and  branch  sequence  numbers 0 

B 0 6  Boundary  Condition  Codification 

A  boundary  condition  of  the  ;nain  member,  or  of  a 

branch  if  any,  consists  of  six  parameters  which,  in  the 

program,  constitute  the  elements  of  a  column  matrix  called 

a  state  vector 0   Five  common  boundary  conditions,  those  of 

a  fixed,  free,  pinned,  propped,  or  a  roller  supported  end 

may  be  indicated  by  a  code  given  below  which  the  program 

can  interpret,  and  from  which  it  constructs  the  required 

state  vector0 
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2 .61   Boundary  Condition  Code 

Codable  boundary  conditions  are  indicated  by  the 
below  listed  values  of  the  floating  point  constant  SVs 
a*   SV=lc   fixed  end 

bo   3V=20   free  end 

i 

Co   3V=3o   pinned  end 

d0   SV=^0   propped  end2 

3 
e0  SV-5o      roller  supported  endy 

«^hen  a  boundary  condition  is  other  than  one  of  the 
five  expressable  by  this  code,  the  user  must  construct  the 
necessary  state  vector  and  cause  it  to  be  read  in  as  data. 
This  situation  will  almost  always  obtain  when  both  Case  1 
and  Case  2  mode  frequencies  are  sought  for  the  same  system0 
It  is  apparent,  for  instance,  that  a  constraint  which  fits 
the  definition  of  a  prop  when  considering  Case  1  vibrations, 
may  not  be  describable  by  any  of  the  codes  available  when 
Case  2  vibrations  are  considered 0 

In  constructing  a  state  vector  it  is  only  necessary 
to  know  whether  the  value  of  a  constituent  parameter  is 
zero  or  non-zero 0   The  code  for  a  parameter  having  a  value 
of  zero  is  0o,  while  that  for  one  having  a  non-zero  value 
is  loo 

The  constitution  of  the  state  vectors  for  Cases  1 
and  2  are  given  below 0 


1  =  3 

•'See  Appo  B,  B064  Catalog  of  Coded  State  Vectors 
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Be62   In-Plane  State  Vector  (Case  1) 

The  vector  elements  in  order  are: 

u  displacement  in  the  x  direction0 

v  displacement  in  the  y  direction 
19*  slope,  in  the  xy  plane,  with  respect  to  the  x  axis, 

M  moment  about  the  z  axis , 

Vy  shear  in  the  y  direction 

N  normal  force  in  the  x  direction 
B063  Out-of-Plane  State  Vector  (Case  2) 

S     twist  about  the  x  axis „ 

T  torque  about  the  x  axis0 

-w  (negative  of)  deflection  In  the  z  direction,, 
"UT  slope,  in  the  xz  plane,  with  respect  to  the  x  axis, 

My  moment  about  the  y  axisc 

Vz  shear  in  the  z  direction 0 


central  axis 


Figo  B063   Vectorial  representation 

of  state  vector  parameters 
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3,6^  Catalog  of  Coded  State  Vectors 
Case  L 


© x 

Boundary  Conditions      Tz 

State 
Vector 
Parameter 

a 

l=r^o 

fixed 

free 

pinned 

propped 

p 

roller  support 

u 

V 

vy 

N 

Oo 
Oo 
Oo 
lo 
lo 
lo 

lo 
lo 
lo 
Oo 
Oo 
Oo 

Oo 
Oo 
lo 
Co 

lo 

lo 

lo 
Oo 
lo 
Oo 
lo 
Oo 

Oo 
lo 
lo 
Oo 
Oo 
lo 

bo   Case  20 


Boundary  Conditions 

%r    x 

State 
Vector 
Parameter 

l_ 

flT^n 

fixed 

free 

pinned 

TTTTTfTT 

propped 

P 

roller  support 

s 

T 
-w 

T 
»y 

Vz 

Oo 
lo 
Oo 

0, 

lo 
lo 

lo 
0„ 
lo 
lo 
Oo 
Oo 

Oo 

lo 
Oo 

lo 
Oo 

lo 

lo 
Oo 
Oo 
lo 

Oo 

lo 

lo 

Oo 
lo 

1. 

0, 

0, 
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B97   Format  of  Data 
Bo 71   Problem  Statement 

The  first  data  card  contains  a  statement  of  the 
problem „   The  card  has  ten  fields  under  the  control  of  the 
following  field  specif ications s 

312  ,  ^F3o0  ,  2F8oO  ,  F2o0 

Field  names  and  the  information  each  field  conveys 

are  given  below  followed  by  an  explanation  of  codes  used. 

Field    Name  and  Content 

number  of  components „ 

number  of  branches 0 

code  calling  for  in  and/or  out-of-plane 
solution 

number  of  modes  sought 0 

code  calling  for  the  use  of  a  lumped  or 
distributed  mass  model  in  analysis 0 

code  calling  for  inclusion  or  neglect  of 
shear  deflection  ef f ects . 

code  calling  for  Inclusion  or  neglect  of 
rotational  inertia  effects . 

upper  limit  of  frequency  range  of  investi- 
gation (radians/sec, ) o 

lower  limit  of  frequency  range  of  investi- 
gation and  initial  interation  frequency e 
If  a  zero  is  placed  in  this  field  the 
search  for  mode  frequencies  will  begin  at 
oOl  radians  per  second 0 

10.   GDABC   code  calling  for  print  out  of  the  frequency 
and  value  of  the  frequency  determinant  for 
each  iteration,  the  number  of  subsections 
assigned  each  component  treated  by  lumping 
mass,  the  boundary  condition  state  vector 
for  both  ends  of  the  main  member  and  of  each 
branch,  the  number  of  components  in  the  sys- 
tem, the  frequency  range  of  investigation, 
the  highest  mode  sought,  and  all  of  the  input 
data0 

B-ll 


1. 

NS 

2. 

NBR 

3» 

I  OP 

^0 

HM 

5» 

~MYK 

6. 

3D 

7* 

HI 

80 

0I4EGAG 

9o 

OMEGAI 

Codes  used  in  the  first  data  card  are  as  follows: 

a.  IOP=0   both  in- and  out-of°plane  mode  fre- 

quencies sought. 

b.  I0P=1   in-plane  frequencies  sought, 

c.  I0P=2   out-of -plane  mode  frequencies  sought, 
do   AMYK=1   lumped  mass  model  desired. 

e.  Ai4YK=0  distributed  mass  model  desired. 

f.  3D=0    consider  shear  deflection. 

g.  SD=1    neglect  shear  deflection. 

h.  RI=0    consider  rotational  inertia. 

i.  RI=1    neglect  rotational  inertia. 

j,  GDii3C=0  print  solution  onlye 

k.  GDaBC=l  print  solution  and  input  data. 

1.   GDa3C=2  print  solution,  input  data,  and 
graph  data. 

B.72   Boundary  Conditions 

Boundary  conditions  for  the  in- and  out-of-plane  case 
are  read  in  separately.   Codes  outlined  in  B.6l  apply  to 
both  cases.   When  both  in- and  out-of-plane  mode  frequencies 
are  sought,  the  boundary  conditions  for  the  in-plane  cases 
are  read  in  first. 

The  field  specification  for  the  first  boundary  condi- 
tion data  card  is: 

22F  3.0 

The  first  N+2  fields  of  this  card  must  be  filled  where 
N  is  the  number  of  branches  in  the  system.  The  coded  repre- 
sentation of  the  main  member's  left  and  right  end  boundary 
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conditions  are  entered  in  the  first  and  last  fields  used 
respectively o   The  codification  of  the  far  end  boundary 
condition  for  each  branch  is  entered  in  the  intermediate 
fields  in  accordance  with  the  branch  sequence  number  (B<,53)° 
When  a  boundary  condition  is  not  describable  by  one  of  the 
codes  given  in  B06l,  a  0o  nmst  be  entered  in  its  data  field. 
For  each  zero  entered  in  the  first  boundary  condition  data 
card  an  additional  data  card  is  required  on  which  is  entered 
the  actual  state  vector  as  given  in  B062  or  B063,  codified 
in  accordance  with  the  last  paragraph  of  Bo6l0   These  cards 
are  placed  immediately  behind  the  first  boundary  condition 
data  card  their  order  dictated  by  the  sequence  number  of 
the  branch  to  which  they  apply  (B053)° 
Bo 73  Branch  Joining  Anprle 

When  the  system  includes  branches,  the  angles  at  which 
the  branches  join  the  main  member  are  entered  in  one  or  two 
data  cards  as  required  whose  field  specification  iss 
20  F  ?03 

For  an  N  branch  system  the  first  N  fields  are  used0 

The  angles,  in  degrees,  are  entered  in  successive  fields 

in  accordance  with  the  "sequence  number  (Be>53)  of  the  branch 

to  which  they  apply, 

B • 7^  Component  Extensive  Properties  and  Branch  Component 
Identification 

A  data  card  Is  required  for  each  system  component. 

The  order  of  these  cards  is  dictated  by  their  component's 


1Field  specification  6F2o0o 
B-13 


sequence  number  (B,51)«   The  card  field  specification  is s 

3F  5,2  ,   2F  6o2  ,   6F  7,2  ,   12 
Field  names  and  their  data  content  are  given,  in  order, 
below s 

D     outside  diameter  (inches). 

2 

DI    wall  thickness  (inches), 

RHO   radius  of  curvature  with  sign  (3,44) 

(inches).  Notes  Applies  to  curved  com- 
ponents only,  A  zero  (0.)  is  entered  in 
all  other  cases, 

TL    length  (inches).   Notes   applies  to 

straight  components  only,   A  zero  (0,)  is 
entered  for  all  other  cases , 

TH2TA  are  central  angle  (degrees )„   Note;   Curved 
components  only0   A  zero  (Cl)  is  entered  in 
all  other  cases 0 

CLX 

CLY, 

CLZ 

CTX 

CTY. 

CTZ 

NImBR  branch  component  identification  number  (B,*t2), 
B,75  Intensive  Properties  of  Components 

The  first  intensive  properties  data  card  has  four 
fields  whose  specifications  are 

F  7o3  ,   2E  10o2  ,   F  2,0 


spring  constants  (3 


m9 


1-2 

These  fields  are  left  blank  in  cards  pertaining  to 

spring  hangers, 

3  These  fields  are  left  blank  except  in  cards  per- 
taining to  spring  hangers. 
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The  field  names  and  their  data  content  are  given  below 0 

AAMU  density  of  first  component  (lbs0/cUo  ft.). 

AE    elastic  modulus  of  first  component  (psi)„ 

AG    shear  modulus  of  first  component  (psi)0 

DD         code,   see  *  below u 

•    When  the  intensive  properties  of  all  components  (spring 

hangers  excluded)  are  identical,  the  DD  field  is  left  blank 

and  the  first  intensive  properties  data  card  becomes  the 

only  intensive  properties  data  card.   If  DD  is  non-zero, 

the  data  on  the  first  card  is  interpreted  as  pertaining  to 

the  first  component  only  and  an  intensive  properties  data 

card  must  be  prepared  for  each  of  the  remaining  components, 

i0e0,  the  second  through  the  last,  spring  hangers  includedo 

These  additional  cards  are  identical  to  the  first  card  in 

format  and  field  content  except  that  the  DD  field  is  deleted 0 

A  blank  card  must  be  inserted  for  each  hanger .   Intensive 

properties  data  cards  are  ordered  in  accordance  with  their 

component's  sequence  number  (B<>51)° 

Bo 8   Data  Deck  Assembly 

The  data  deck  is  ordered  as  follows: 

lo   problem  statement  card, 

2,   boundary  condition  data  card(s)0 

3«   branch  joining  angle  data  card(s)0   (deleted  when 
system  has  no  branches )0 

*K   component  extensive  properties  data  card(s)0 

5°   component  intensive  properties  data  card(s)0 
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The  user  may  obtain  solutions  for  any  number  of  prob- 
lems desired  by  "stacking"  assembled  problem  data  card 
decks  at  the  end  of  the  program  deck.   Note  that  a  blank 
card  must  be  inserted  between  the  last  card  of  the  program 
deck  and  the  problem  statement  card  of  the  first  problem 0 
A  blank  card  must  also  be  placed  behind  the  last  intensive 
properties  data  card  of  the  last  problem,   DO  KOT  insert 
blank  cards  between  problems.   The  problem  statement  card 
of  the  N+first  problem  must  be  immediately  behind  the  last 
intensive  properties  data  card  of  the  IvJtn  problem. 
3.9   Flexibilities  of  the  Program  VIPIPE 
B.91   Shear  Area  Form  Factor 

A  shear  area  form  factor  of  one  half  is  built  into 
the  program.  L&j   It  may  be  altered  by  changing  card  OOO58Q0 
B.92   End  Conditions 

A  properly  constituted  boundary  condition  state  vector 
must  have  three  zero  and  three  non-zero  elements.   Fre- 
quently, however,  an  actual  end  condition  may  not  meet 
this  specification,  e.g.,  an  end  of  the  main  member  or  of 
a  branch  may  be  supported,  or  connected  to  a  piece  of. 
machinery  in  such  a  manner  as  to  approximate  being  held 
in  a  rubber  mounting  block  (Fig.  B.92-1),   To  program  this 

fictitious  free  end -^> 

system-^ 

Fig.  92-1   Fictitious  end  condition 
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end  condition,  the  user  should  make  a  good  approximation 
of  the  elastic  properties  of  the  mounting,  support,  or 

connection  in  terms  of  the  linear  and  torsional  spring 

2 

rates  of  the  program  (Bc^5,  Bo7^),  then  attach  a  fictitious, 

short  (perhaps  ,01  in„)  massless  straight  section  to  the 
end  under  consideration,  and  specify  in  the  boundary  con- 
dition for  that  end  that  the  end  is  free0 

Through  the  use  of  fictions  of  this  sort,  properly 
constituted  end  condition  state  vectors  may  be  constructed 
for  any  real  situation  with  little  or  no  loss  of  accuracy,, 
B.  93  Subsectionlnp; 

Subsectioning  for  the  lumped  mass  treatment  of  com- 
ponents is  carried  on  in  subroutine  SUBSEC  «   The  subroutine 
is  set  up  to  divide  a  component  into  a  maximum  of  ten  sub- 
sections.  The  writer,  through  the  analyses  of  several 
systems,  found  that  this  figure  gave  good  accuracy  of 
solution  for  a  reasonable  expenditure  of  machine  time0 

For  all  components  having  a  length-to-dlameter  ratio 
greater  than  six,  the  number  of  subsections  is  computed 
as  two  times  the  value  of  a  parameter  HM  which  appears 
in  the  subroutine 8   The  user  may  control  its  value  by 
altering  card  00^720  in  the  subroutine » 


1See  Appo  A,  sections  A022,  and  A032« 
2 i.e.,  treat  the  mount  as  a  hanger . 
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B.9^   Values  of  the  Frequency  Determinant  and  Solution 
acceptability  Criterion 

The  frequencies  for  which  the  value  of  the  frequency 
determinant  is  zero  are  the  natural  frequencies  of  the 
system  which  it  describes c   In  light  of  the  number  of  sig- 
nificant figures  available  in  a  computer  and  of  the  effects 
of  round-off,  however,  insistance  that  iteration  and  con- 
vergence operations  continue  until  the  value  of  the 
frequency  determinant  actually  became  zero  would  require 
an  excessive  expenditure  of  machine  time.   Instead,  itera- 
tion proceeds  until  the  sign  of  the  value  of  the  frequency 
determinant  changes,  then  a  convergence  routine  is  carried 
out  until  the  difference  between  the  values  of  frequency, 
for  which  the  corresponding  value  of  the  frequency  determinant 
have  opposite  signs,  is  less  than  a  specified  value  CR, 
The  result  of  a  linear  interpolation  between  those  fre- 
quencies is  accepted  as  a  natural  frequency  of  the  system. 
The  specified  frequency  difference  CR  is  termed  the  solution 
acceptability  criterion,, 

To  compute  a  starting  frequency  increment,  the  program 
constructs  a  synthetic  straight  pipe,  equal  in  length  to 
the  length  of  the  main  member,  and  having  diameter,  wall 
thickness,  and  intensive  properties  equal  in  magnitude  to 
a  weighted  average  of  those  properties  of  the  composite  main 
member.   The  fundamental  frequency  of  the  synthetic  pipe 
is  computed,  and  one  half  of  its  value  is  taken  as  the 
startinp-  frequency  increment 0   The  initial  solution 


This  increment  may  be  changed  by  altering  card  001090* 
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acceptability  criterion  is  taken  as  one-twenty-fifth  of 
the  starting  frequency  increment e  The  initial  solution 
acceptability  criterion  may  be  changed  by  altering  card  001110. 

The  starting  frequency  increment  and  the  initial  solu- 
tion acceptability  criterion  are  used  until  the  second 
natural  frequency  is  found,  after  which  the  frequency 
increment  is  recomputed,  as  it  is  after  each  succeeding 
natural  frequency  is  determined,  and  given  a  value  equal 
to  one  half  of  the  difference  between  the  last  two  fre- 
quencies found o   The  new  solution  acceptability  criterion 
is  then  recomputed  as  one-fiftieth  of  each  new  frequency 
increment .   The  new  solution  acceptability  criterion  may 
be  changed  by  altering  card  002800 « 
B.95  Flanges  and  Couplings 

Pipe  couplings  and  flanges  represent  real  "lumps"  of 
mass  mounted  in  a  piping  system 0   They  may  be  treated  as 
individual  system  components.   If  this  is  desired,  a  com- 
ponent sequence  number  and  a  branch  component  identification 
number  should  be  assigned  in  the  usual  manner  (B.51.  and 
B.52),  and  appropriate  extensive  and  intensive  properties 
data  cards  (Bo?**,  and  Bo75)  should  be  included  in  the  data 
decko  These  components  are  treated  in  the  program,  as 
rigid  sections  having  distributed  mass. 
B„96  Hangers  Exhibiting;  Coupling 

Hangers  such  as  that  shown  in  Fig0  B.96-1  can  be 
treated  with  the  tools  available  in  VIPIPE  by  considering 
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the  mounting  arm  as  a  branch  having  a  hanger 


Fig,  B.96-I   Cantilevered  hanger 
and  a  free  end.   Recourse  to  similar  methods  will  permit 
treatment  of  almost  any  hanger  configuration. 
B.97  Systems  Having  Unusual  Flexibility 

When  the  main  member  of  a  system  is  essentially 
cantilevered  the  computed  starting  frequency  increment  may 
be  so  large  that  the  first  two  system  natural  frequencies 
are  not  found.   For  analysis  of  such  systems  it  is  recom- 
mended that  the  constant  ^.7300^019  in  card  001080  be 
changed  to  I.875IOO67,   Similarly,  if  it  is  suspected  that 
the  fundamental  frequency  of  the  system  is  less  than 
•01  rad./sec.  a  suitably  lower  starting  frequency  should 
be  programed  (B.71). 
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APPENDIX  C 


FLOri  DIAGRAMS 


Col  General  Remarks 

Flow  diagrams  are  included  for  all  parts  of  the  main 
body,  and  for  several  subroutines „  Those  subroutines  for 
which  no  flow  charts  are  included,  with  the  exception  of 
INVERT,  incorporate,  at  most,  two  decisions  in  their 
structure*   as  applicable  they  sense  on  SD  and/or  RI  (  A»31) 
to  change  the  value  of  one  or  two  elements  of  the  transfer 
matrix  constructed  by  the  particular  subroutine e 

Subroutine  INVERT  performs  a  matrix  inversion  by  a 
pivotal  technique .   Sufficient  comment  statements  are 
included  in  the  subroutine  to  define  its  structure 0   If  the 
user's  FORTHaN  library  includes  a  satisfactory  single 
precision  matrix  inversion  subroutine,  INVERT  may  be  deleted 
after  suitably  modifying  the  dimension  (cards  010070, 
010080  BRANCH,  and  011080,  011090  BRaNCO)  and  call  (card 
010230  BRANCH,  and  011660  BRANCO)  statements  of  subroutines 
BRANCH  and  BRANCO . 
C«2   Index  of  Flow  Diagrams 
C.21  Main  Body 

Flow  Diagram  Page 

INPUT  C-3 

INVARIANTS  C-5 

CONTROL  IN-PLANE    1  C-6 

CONTROL  OUT-OF-PLANE  C-6 

ITERATION  C-7 

OUTPUT  .           C- 9 

1See  C03  c-1 


C.22  Subroutines 

Flow  Diagram 

SUBSEC  C-10 

MATMUL  C-  8 

BRANCH  C-ll 

BRANCO  C-12 

STAVEC  C-13 

STAVEO  C-13 

FINMAT  C-l^ 

FINBRa  C-l^ 

DELMA  C-14 

C.3   Flow  Diagrams1  c~3  through  c.^ 

Note:   Since  CONTROL  IN-PLANE  and  CONTROL  OUT-QF-PLANE 
are  identical  in  structure,  only  the  former  flow  diagram  is 
included.   They  differ  only  in  the  subroutines  called  and  the 
state  vectors  used  (A. 22).   Subroutines  SUBSEC,  MATMUL,  INVERT, 
FINMAT,  FINBRa,  and  DELMA  are  common  to  both  cases. 


*A  complete  list  of  symbols  is  given  in  Appendix  A, 
section  A. 3. 
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\hEAD:    Dimension? 
\stat,ement  / 

. JL 

\READ:    Problem/ 
Statement / 


\READ:    Branch   joining? 
Wgle     PHI(I) 

I  =  1,NBR 


p 

(l),DI(l),RHO(I), 


READ:    D(I),DI(I),RHO(l),TL(I),THETA(I)v 
\CLX(I),CLY(I),CLZ(I),CTX(I),CTY(I), 
\CTZ(I),NNBR(I) 

1=  l.NS 


\READ:    AMU(I),E(I),0(I)/ 
\ I=2.N8  / 


J'.'PT'T_ .  '-low  .  Tju  y-n\ CO\; 


^ 


AMYK«AWYK-1. 

OME'iAO  =OME^AI 

FFAC  =  .  5 

L  =  l 

M  =  l 

AA-  PR  =  CC  =  OC=-EE=0. 

AaMTT=AL=AD=AID=0. 

KKsO 

NNN^NPR-Hl 


THETR(I)=*THETA(I) 
DT(I)  =  DI(I) 

RAMU(I)  =AMU(I) 
DI(T)  =  DT(I)    -   2.0-m-DKI) 
1=  1,NS 


|    GO  TO   I N VARIA NTS 


_INVfeHIAl\TS      FLOW      DIA  >HA?1 


>0       X  ^ix(i)frLy(i)*-CLZ(i)-P 

— <      C'PX(  I  )  +  CTY(T)+-CTZ(  I) 

\  1  =  1, NS 


|AMU(I)=AMuTl)»Tr  #(D[f J**2-DI ( I )**2 )/( 6912 .#32 .  17*18 .  )] 


kHETA(l)  -  THETAU  )*ir/l80^       1  ^  <V, 

|tl(i  )  =  -n-;sTA(  i  )»afsf(fho( i ) )  I 


AJT( I )=TT  #(D( I )»*4-DI( I )*#4 )/38. 
AJ(I)=  A.rT(I)/2. 
P(I)-1./(A.T(I)#E(I)) 

AIX(I)=  SQRTF( (D(T)**2    HI ( I )*#2 )/8. ) 
AIV(T)=AIX(_I)/1. 414214 


W 


Ti*P  AJhTANTO 


<&> 


AL=  AL-t-TL(T) 

AD  =  ADfb(I)»'T'L<T) 

•  E  =  EE+E(I)*TL(I) 
AID  =AID+DI(I)*TL(I) 
AAiW  =  AAM1H-AM1T(  1  )*i  L(  I ) 


-4SH> 


AD  =  AD/AL 

AID  =  AID/AL 

AFR=  1T/64 . # ( AD**4-AID#*4  ) 

FR  =  4,7S*#2*SQRTF(EE*AFR/(AAMU#AL##4 ) 

AD0M=PR/2.0 

DOM  =  ADOM 

CR  =  PR/25, 

PDOM  =ADOM 


[TO  CONTROL   OUT  OF  PLfvNC~| 


PI,A':E  FLOW  DIAGRAM 


JFr  "J  A   INV, 
\ 


-{CALL   STAVECi 


ISVI(J,1)  --.  SVIP(J,l)h 

I  '  .£-1.6  I 


PSVE(  J.l  )  =  SVIP(  j",NMPH  f^*— -<'SV(NMFF  )  >-' I  CALL  '"STAVECi 

""•""[jk-o}' 

Sit] 

<SV  (N  )J>~^         *-{CALI,   STAVECi 


|CALL  RIGID |  |CA~LL  DISTM]  fCALL   STIFCO"[-^-<Z(I :)>^-fCALL   CFIELD]-*-pHCALL  SFIELD>«Az(  I \ 


| CALL   POINT  1 
[CALL  MATMUlI 


I^ALL  BFA1.CH  \~ 


[to  iteration! 


1 


_ITHHA'?'TnN FT.OW  LIA''.:-AM 

|  CALL  DELMA] 


FJDtL=  DEL(L) 
OMEGA  (L+D-OME; 
<0 


I     [Numerical   val 
[determinant   c 
,HA(L)-t-2..-*nR  I 


1  TO  OUTPUT] 


OMEGA C  =  OMEGA  (L-l) 

AA=  1. 

BB  =  1. 

DELC  =  DEL(L-1) 


-^2-<^MEGA(L)-0MEGA(L   1)-CR> 
>0 


ADOM=  (OMEGAM(M-l )-0MEGAM(M-2  )  )*.( 
CR=  ADOM/50. 


|  DOM  =  ADOM 
~1 


?  output] 


3 


|  OlviESha»H  =  OME0.AM(M-l\  +  2.»CrT| 

jE 


fe^^^f]  (5>H  omegao^omegaVl-h)  | 

[goto   AJ<J,esri  loozl      I 
I  Com  ikoi    car-aP-PLftHE  | 


(l)-,9*dom| 

. 5*DOM 


S'TBHO'tTIN  E      KATM'IL      RT.OW ni 


r  cu,K)=i;   I 

J  =  l,6     K-1,6 


|D(J,K)  =  D  (  J  ,  K  )+A(  J ,  L  )*U  ( L ,  KJl 
L        L  =  1,6 K  slt6     J  =  1,6     j 

, I , 

IT(.T,K)  =D(J,K) 
[,1  =  1,6     K  =  l,6  j 


OUTPUT      FLOW      lJ±<y_LkM_ 


P=  JJ4-I] 


p(T)^o. 


MN  =  AA  =  BB  =  CC  =  GG  =  KK  =  0. 

ADOM  -  BDOM 

CR  =  BR 

OMEGA ( 1 )  =  OMEGA I 

OMEGA 0  ■ OMEGA I 

DOM  =  ADOM 


GO    TO  Address    450  "of] 

[input 


CONTROL   OUT-OF-PLANE 


U   E  .  EkAiLGH     FLOW 


iifFl 


i  ')■ 


I  Construct    third   order   coordinate    j 
I  transformation   matrides    "rl    ar-d   02j 

I  RlIJ,K)=  0.  I 

Tsl.3      K=-  lj_3_J 

ZI1 , 

1(J,K)»  R]  (J,K)+^(.T,L)#HfL,Kn 
L  =1,3      K ■_=  1 ,3      J  =  1 13  J 


_     -         I-         , 


pi(i,i )  =  n(?,"n )  =  '"3,3  j  =  1.1 

|n(4J4)=  M(5,b)=  7T_(^,6|=  1.  j 


, I ... 

!  TI(5,L)=  R(2,L) 
j  I'(6,TJ  =  R(3,L) 
!  L  ■  1,3 


I       R(.T,K)=0.       1 

j  T  =1  ,5     K  -  1 ,  3  | 


S;rPRO:jTT' E.  _F'FA!'.CO   FT  OV;  pik'-ll-X 


. -■    -i.~_ -    

Reorder  boui.dt.ry  condition  state 

vector  such  that  the  first  three 

elements  pertain  to  displacements, 

I  the  last  three  to.  forces. 


p-eorder  branch  transfer  matrix"! 
I  to  make  It  compatible  with  the  j 
I  reordered  state.,  vector,. 


r .1    \ 

L.T=1^3   K=l,3 

r..      .1         ■ 

rS(J,K)  =  S(J,K)+Rl(J,L)*Gl(L,Kl] 
; L=  1,3  K  =  1,3  Jg  1,3 


U(l,l)=  U(2,2)=U(3,3)=  1 
U(4,4)=  U(b,5)  s  U(6,6)   " 


U(2,l)=  S(l,l) 
U(2,3)=  S(l,2) 
U(2,4)-  S(l,3) 
U(5,l)=  S(2,l) 
U(5,3)=  S(2,2) 


1. 

U(b,4)=  S(2,3) 
U(6,l)  =  S(3,l) 
U(6,3)=  S(3,2) 
U(6,4)=-  S(3y3) 


C~\Z 


S'TBROT'T'KE      S'l'AV/r      '-'LOW      DTA   ^AM  . 


{_RC_(4,N  )  =  f.C(fa,N  )  =  BC(6.N  )  =  1.  - 


/av-p>>sQ 


r 


H 


5V-5^>— Hpr(^,"w  >  =  M [5 j,N  )  a  PC ( 6 ,H  ) £l  [1- 


=  0 


-^return] 


■»-{pC~("i,N)  =  PC('a,N~)  -BP  (  3  >  )_  =•  1 7}— H"-E^ 


'■P-N'1 


<SV-4.>-^- *-|;PCU/'J=  r>'r  I  z\  r\  =s  PC]  5AN  )  =T7] HRBTitpm  ] 

-->-{p72,/\- )  =  P^  ( 3,_M  j  =  PC  ( 6,n~[=  lV}- HZE^i'] 


S-PRQTJ'PINE_  STAVEO     PLOW 


Sl'V-r-.CVlTJNE      Fi    !v>.T      i-LOVn      MA  :>.A)V. 


• ;  I- 1>  =  °         *-fv  (  J,  K")  =  U(  J,  K)'\  .-—-*-)  RETURN  ] 

N^  |    J-  1,6  K  =  1,6  I 


V(J.K)=  0.        i- 
.7^1,6     K-1,6  1 


a-    V(.I,K)  -  V(  J,K)+U(J,L)*UU(L,K)  j- 
L=  1,6      K  =  l,6      J-  1,6  j 


iUBROUTDJE     PINFRA      FLOW      DIAORAM 

—^return! 


V(J,K)«U(J,K) 

<TINFR-1  V^— •JjT  =  1,6      K  =  1,6 


V(  j7k)=  0.    ""V— *-T\'(J,K)="*V(J,K)+U(  J,"l)*VV(L,K)  h- 

r»l,6     K  aj.6  I  L  _  I»"l»ft K  =  1.6     J=1.6 j 


— *- Treturn 


I     M~0| 


_::r_ 


i 


j-6>-^^-H  j=  j* it- 


^— <SVE(N,1  J   >»- 

sect 

p7TL,K)=  vv(n7kT1 
L k»i.s       I 


aPi'ENDIX  D 
PROGRAM  LISTING 


D.l 


Index 


Dell   Main  Body 

Section 

Sequence  Numbers 

Page 

INPUT 

000040-000760 

D-2 

INVARIANTS 

000790=001170 

D-3 

CONTROL  IN-PLANE 

001200=001810 

D-4 

CONTROL  OUT-OF-PLANE 

001840-002450 

D=5 

ITERATION 

002470-003260 

D=6 

OUTPUT 

003280-004660 

D-8 

Do  12   SUBROUTINES 

SU3SEC 

004690-005140 

D=10 

DISTM 

005170=005710 

D-ll 

DISTMO 

005740=006280 

D-12 

SFIELD 

006310=006640 

D-13 

SFIELO 

006670=006910 

D=14 

CFIELD 

006940=007370 

D=14 

CFIELO 

007400-007890 

D-15 

POINT 

007920=008090 

D-16 

POINO 

008120-008310 

D-17 

MATMUL 

008340=008900 

D-17 

RIGID 

008670=008643 

D-18 

RIDIO 

008930=009160 

D-18 

STIFCO 

009190-009410 

D-18 

STIFOO 

009440-009660 

D-19 

HANGER 

009700=009850 

D-19 

HANGEO 

009880=010030 

D-20 

BRANCH 

010060-010730 

D-20 

STAVEC 

010760=011040 

D=21 

BRANCO 

011070-012210 

D-22 

STAVEO 

012240-012480 

D-24 

INVERT 

012510=013220 

D-25 

FINMAT 

013250-013380 

D=26 

FINBRA 

013410=013540 

D=26 

DELMA 

013570=013750 

D-27 

D,2 

Listing 

D-2  througl 

1  D-27 

D-l 


PPCCR/1"  VIPIHl 

E  0C0C10 

C  CCC030 

C      INPUT  OOCCUC 

L  OCCC5C 

C  CCCC60 

NM«0  0CCC70 

JJ=i  OOOCRC 

M500DIMLNSION  4JI  (5C),4JI  50),».(  :C)  ,MX(SO)  ,AIY(  50)  ,PHII2C)  .SVI  (b,  1),     0CCC9C 

1SVB{6),D(50),CI<5C),TL<5C),Z(50),G(5  0),E(50),AMI(50),THETA(S0),RH0  .000100 

2(50)  .CN-tCAV^O)  ,  C  EL  (  SCO  >  ,CM  t  CA  (  5CC  ),  A(6,6),P{6,6),U(t-,6)  ,CLX(5C) ,C   0C011C 

3LY{5C),CLZ(5G),CTX<50) ,CTY( 50) ,C1Z(50> , VI 6, 6 ) , VV( 6,6) ,UU (6,6),SVt(   0  001 20 

«*6,1)  ,SV(2C),SVIP(*,22),SVP(22),SvnPC  6,2>  ).NNBR(50'),6AI"U(-6C)  ,THETR(   0C013C 

520) ,DT(5C),P(20)  000140 

REAL  15CC,NS,NeR,I0P,HN,AyYK,SC,'U,,DMECAG,0M:CAIfG0ABC  0C01 bC 

1500  FCRMAT(  3I2,4F3.C,2F8.C,F2.C)  000I6C 
NyBR=NBS*2  0C017C 
IF(  ICP-1 ) 1501,1501, 15C6  0C0180 

1501  REAC  1502,  (  bV  (I)  ,  1=1,  NKPP)  000190 

1502  FCRfAT(22F3.0)  00020C 
CC  1505  J=1,NNBR  0C0210 
IF(SV( J) ) 1503, lr03, 1505  000220 

1503  REAL  15C4,(SVIP<  I,  J  ),  1  =  1  ,6)  OC0230 

1504  FCR^AT(6F2.0)  000Z4C 

1505  CONTINUE-  0C0250 
IF(  IOP) 15C6, 150C, 1510  000260 

1506  READ  15G7,(SV'C(I)',r  =  l,M«ER)  000270 

1507  FCRMAT(22F3.0)  0CC280 
CC  1510  J=l,NPfi3  0C0290 
IF(SVCU)  )  1508,  1508,  1510  0C0300 

1508  REAL  1509  ( SVCP (  I , J  ) ,  I  =  1  ,  6  )  000310 
15C9  FCRMAK6F2.0)  000320 
1t>10  CONTINUE  000330 

IF(NER) 15  13, 1513, 1511  OC0340 

1511  RFAO  1512, (PH1I  I  ),  1  =  1, NBR)  000350 

1512  FGRVAT  (1CF7.3)  000360 
15130READ  1514,(01  I  )  ,  C  I  (  I),RHC(  I  ),TL(  I  )  ,T HCTA ( I ) ,CLX ( I ) ,CLY(I  )  ,CLZ(I),C   0C0370 

1TX( I ),CTY(  I  ) ,CTZ( I ) ,NNPR (I  )  ,  1=1  ,f\S)  000 38C 

151U  FCR^AT(iF5.2,2F6.2,6F7.2,I2)  00039C 

REAU  1515  ,  AAPU,AE,AG,CC  O0C4CC 

1515  FORMAT  (F7.3,2E10.2,F2.C  )  .  000410 
IF  (  DC ) 15 16 , 1 5 16 , 151 e  000420 

1516  DC  1517  1=1, NS  000430 
AMU(  I  )=AAN"U  000440 
E(I)=AL  O0C450 

1517  G(  I)=AG  000460 
GC  IC  1520  00047C 

1518  AMU(  1  )  =  AAMU  0C048C 
E( 1  )=AE  0C0490 
G(  1  )=AG  000500 
READ  1519, ( AMU( I ),EU ) ,G( I ) , 1  =  2, \S)  00051C 

1519  FCRVAT  (F7. 3.2E10.2)  0C0520 
IFICPCGAI ) 1520, 1520, 1521  000530 

1520  GVEGAI=.01  000540 

1521  C^EGA( 1 )=CMtGAI  00055C 


AMYK=AyYK-l.G  OCCitC 

CMEGAC=CMEGAI  OC057C 

FrAC=.c  0C058C 

L  =  1  0CC^90 

y=1  0CC6CC 

AA=C.O  00061C 

BF=0.0  0CC6'0 

CC  =  0.0  OC06-50 

GG=0.0  OCCfcUO 

KK=C  00G65C 

AAMU=G.0  0CC660 

EE=0.0                                                    -  000670 

AL=0.0  OC0660 

AD=O.C  0C069C 

AID  =  C0  0C070C 

NNN=NPR+1  0C0710 

DG  1522  1=1, NS  OC072C 

DT  (  I  >  =i: I  (  I  >  0C073C 

BAMUm  =  AI*U(I  )  000740 

THE1RI I )=THETA( I )  0C075C 

1522  Did  )=C(  I  )-2.C»CH  I  )  00076C 

;  0C0770 

;  000730 

Z                  INVARIANTS  0C0790 

:                                       ,  ocoecc 

DC  1C8  1=1, ns  ocoeio 

IF(CLX{l)+CLY(l)+CLZ<I>+CTX(l)+CTYU)+CTZ(I))108,101,1C8  0  0C82C 

101  AMU(  I)=AMU(  I)«3.1415926S4*([:m»«2-DI(  I  )  *»2  )  /  (  69  1  2  .♦  32.1  7»  1  2  .  )  0008  30 
IF(RHC( I  )  )  102, 103, 102  000840 

102  THE1  A(  I1=THlTA<  I  )«3.1415<;2654/180.  000850 
TL(  I  l=THETA(I  )»ABSF(RHC(  I  ) )  000860 

103  AJT( I  )=3.141592654*(0( I ) «*4-CI ( I )»»4  )/3I.  000870 
AJ( I  )=AJT( I )/2.  000880 
R{ I )=1 .0/(AJ( I )«t( I ) )  00089C 
AIX(I)=SGRTF( (C(I)«»2+DI< I)«»2)/P. )  000900 
AIY(I)=AIX(  IJ/1.U1U21U  000910 
IF(AMYK) 104, 1C5, 1CU  000920 

104  Z(I)=1.0  000930 
IFJRHCU I ) 1105,106,105  0C09U0 

105  HHM=HV  00095C 
CALL  SUSSEC  (RHClD.THEraj  I),D(I),TL  (D.Hf.ZUn  000960 
HW=HHN  000970 

106  IFINNPRI  I  )  ) 1C8, 107, 108  000980 

107  AL=AL+TL(I)  000990 
AD=AC+H(I)*TL( I )  001COO 
EE*EE+EII1*TL{I1  001010 
AID=AIC+DI( I )«TL< I)  001020 
AAML=AANU+AMU( I )«TL( I )  001030 

108  CCNTINUE  001C40 
AD=AC/AL  0C1C5C 
AIO=AIC/AL  001060 
AFR  =  3.1U159  265l4/6U.«(AC»«U-^m«*i,  )  0010  7C 
FR=4.730040  19««2»SURTF(EE*ArR/(AAMU»AL»»4))  0C1C80 
AL)0V  =  FR/2.0  001090 


DOK=ADCM 
CR=FR/25.C 
BDOI«  =  ACCM 
BR=-LR 

or:  ico  m  =  i, npr 

109  PHI  {  I  I  )=PhI  ( I  I )»3. 141 592654/ 180.0 
IF(  IGP-1  )  TIC, 990, 998 

1 10  MN=1 


L      CONTROL   IN  PLANE 
C 

999  IF(SV(1) 1101 1, 1C21,  1G1  1 
1C11  CALL  STAVLC  (SVI1I,  1.SV1  ) 

GO  TC  1GU1 
1021  DC  1031  J=1,6 
1031  SVKJ.l  )=SVIP(J,  1) 
1CU1  IF(SV(NNPP.)  )  1051  ,  1061,  1C5  1 
1C51  CALL  STAVECt SV(NMPR),  l.SVC) 

GO  TC  1081 
1061  UP  1G71  J=l,6 
1071  SVE( J, 1 )=SVIP( J.NMPR) 
1081  JK  =  C 

IF (NTR) 1000,1 CCC, 1085 

1085  CO  1087  N=2,NNN  , 
IFlbV(N)  ) 1087, 1087,  10^6 

1086  CALL  STAVEC(bV(N),N,SVIP) 

1087  CON  I  INUh 
1000  N=2 

11=1 

DO  10C1  1=1, NS 

P(I)=CLX<  I)+CLY{  I)+CTZ(  I  ) 

IF(P( I) 137,38,37 

37  CALL  HANGER(CLX( I),CLY(  I  ),C1Z(  I)  ,U) 
GO  TO  62 

38  IF(AMYK)39,41 ,39 

39  1F(KHC(I ) >40,42,4C 

40  IF(Z(I))43,45,4? 

41  IF(RHC( I  )  140,48, 4C 

42  IF(Z(  I  1)46,47,46 

43  CALL  CFIELD  ( R ( I ) , Z ( I ) ,G (I  )  , SO ,RHC ( I  ) , THET A (  I  ) , D (  I ) , D I ( I  ) , E < I ) , A ) 

44  CALL  POINT  ( AM U (  I  1 , AI Y ( I  )  ,  O EGA ( L ) ,T L ( I )  ,Z (  I  1  ,  R  I  , R  ) 
CALL  MATMUL  (  A,«,Z(  I.J.U1 

GC  TC  62 

45  CALL  ST1FC0  { THETA(  I  )  ,RHC(  I  ),U) 
GO  TO  ft2 

U60CALL  DISTV  ( AMU ( 1 ) , A  I Y (  I  )  ,GMEGA ( L ) ,T L i   I ) ,R ( I ) , 0 (I) , D  I  (  II  , G (  I  )  ,  SC,R 
1I,E(I),FFAC,U) 
GO  IC  62 

47  CALL  RIGIC  { AMU (  I  ) , TL ( I ) ,CM EGA ( L ) , AI Y ( I )  ,RI , L ) 
GC  TL  62 

48  IF(ZU)  >49,47,49 

49  CALL  SFIELD  ( U I (  I  )  , TL (  I  )  ,  R (  1  1 , Z (  I  )  ,3 (  I  1  , SD, D (  I  ) , E ( I  )  , FFA C , A ) 
GC  TC  44 


0C1 ICC 
0C11  10 
0C1  12C 
001 1 3C 
001  14  0 
0C1 150 
001  160 

ooi  i re 

0C1  lrC 
0C1  IOC 
0C12CC 
001 2  1  C 
001220 
0C1230 
0C12«C 
00125C 
001260 
001270 
0C128C 
00129C 
OC1300 
001310 

oci^rc 

00133C 
0C134C 
OC1350 
0C1360 
00137C 
00138C 
001390 
0C14CC 
001410 
001420 
001430 
00144C 
001450 
001460 
001470 
00146C 
001490 
0C1500 
001 'j  10 
00152C 
001530 
0C1540 
00155C 
001560 
001570 
001580 
0C1590 
001600 
0C1610 
0C162C 
001630 


62 

lF(NNPR(I))65,6-<,65 

63 

CALL  FIN*'AT  (L,IU,I,V) 

DC  c4  J=l ,6 

HO  64  K=l  ,6 

64 

UU(  J,K)=V<  J,K1 

gc  rc  1001 

65 

CALL  FIN?RA  (  L,  VV,NNr.R(  I  )  ,  V  ) 

DO  66  J=l,6 

DO  66  K=l,6 

66 

VV( J,K)=V( J,K) 

IF(NN3R(  I  )  W8,  67,67 

67 

IF(NNbR{  I  >-:>)  1001  ,68,66 

68 

CALL  6RANCH(SVIP,VV,PHI( II > , 

N=N+1 

11=11+1 

cr  tc  63 

1001 

CONTINUt 

GO  TC  2000 

CONTROL   CUT  CF  PLANE 

998 

IF(SVC(1 ) 11013, 1023,1013 

1013 

CALL  STAVEO(SVC( 1  ),  l.SVI  > 

0C1640 
0C16r.0 
00166C 
0C1670 
001680 
0C1690 
001700 
00  17  10 
0C1720 
001730 
001740 
001750 
001760 
00177C 
00178C 
OC1790 
0C18C0 
001810 
001820 
0C1830 
001840 

ooiesc 

001360 
0C1870 

go  to  1043    '  '  ooieec 

1023  00  1033  J=l  ,6  001890 

1033  SVI ( J,l )=SVCP( J, 1 >  001900 

1043  IF(i>VO(N^BR)  )  1053,  1C63,  1C53  CC1910 

1053  CALL  STAVEOISVO(NNBR)  , 1 ,SVE)  001920 

GC  rc  1083  i  001930 

1063  DO  1073  J=l,6  001940 

1073  SVE(JtU=SVCPCJ,NP6R)  001950 

1083  JK=1  001960 

IF(NF.R) 1002, 1C02, 1093  001970 

1093  DO  1095  N=2,NKN  001990 
IF(SVO(N111095, 1095, 1094  001990 

1094  CALL  STAVEC(SVC(N) ,N,SVCF)  002000 

1095  CONTINUE  002010 
1002  N  =  2  O02C20 

11=1  002030 

DO  1003  1=1, MS  002040 

P(  I  )=CTX(I)+CLZ( I)  +  CTY< I )  002050 

IF|P(III1,2,1  0C2060 

1  CALL  HANGEO<CTX(  IJ,CLZ( I  ),CTY{ I) ,U)  002070 
GC  TC  14  002C80 

2  IF(AKYK)3,11,3  0C2090 

3  IF(RFC(I))4,3,4  002100 

4  IF(Z(I))6,5,6  002110 

5  CALL  STIFCO  ( THF T A {  I  ) , RHT (  I  )  , U  1  002120 
GC  1C  14  002H0 

60CALL  CFIELO( AJT(  I) ,R(  I  ),Z(  I  ),G( I ) , S3,RRC(I 1 , THET A (  I  ) ,D( I  )  ,01 (  I  ) , FF   002140 

1AC.A)  OC2150 

7  CALL  PCINC  (A^U(  I), AIX(  I  ),AIY(  I),CMEGA(L),TL(1),Z(  I)  ,RI,B)  002160 


Gr  rr  l^ 

0C21>?C 

8 

I  F  (  L  (  I  )  )  1  0  ,  9  ,  1  0 

0C219C 

9 

CALL  RICIC  (AfU( I),TL( I ) ,AIX{ I ),CMEGA(L) ,AIY(I ),AJT( 1) ,R I,U) 

D02..CC 

go  rr.  iu 

0C2_  IC 

100CALL  DISTNOJ  AMUUlf  AIX(I).f  Al.Vl  t)  ?'OMEGA{L)-fTLiM  i  AJT(I)  #1  (I)  tp(I)  fD 

CC222C 

' 

1  I(  I  )  ,G(  I  )  ,SC,RI,FFAC,n> 

OC2*3C 

GC  IC  Hi 

0C2^UC 

11 

IFUI-C<  I  )  )U  ,  12, U 

OC2250 

12 

I  F  (  Z  (  1  )  )  1  3 ,  9  ,  1  3 

002^60 

13 

CALL  SFICLO  [CI(I),TL(I) fAJT( I»,R( I) ,Z(I 1 tG( II,SC,D( I) ,FFACtA) 

00227C 

Gr  TL  7 

0C22TC 

1U 

IF(NNhR(  I))  17,lr;,  17 

0C2^90 

15 

CALL  FINNAT  (U.UL,  I  ,V) 

0023C0 

DO  16  J=l,6 

0C231C 

DO  16  K=l  ,6 

OC232C 

16 

ULM  J,K)=V( J,K  ) 

OC2330 

GO  IC  1C03 

002340 

17 

CALL  FINBRA  ( L , VV , NNPR ( 1  )  , V  ) 

0C235C 

DO  18  J=l  ,6 

002360 

DO  18  K=l,6 

0C237C 

18 

VV( J,K)=V(J,K) 

0C238C 

IF(NNRR( 1)132,19, 19 

00239C 

19 

IF(NNPR(  I  )-3)  1CC3,82,«2 

0C2U00 

82 

CALL  BRANCO(SVCP,VV,PHI( 11 ) , N , 0 ) 

0C2M0 

N=N+1 

002420 

11=11+1 

0021430 

CC  IC  15 

002iii*0 

1003 

CONTINUE 

0C2H50 
OC2146C 

ITERATION 

002U7C 
0C2UQC 

2C00 

CALL  DELPA(UUf SVI.SVE, V) 

0021490 

0DEL(L>=V(  1,  1  »«V(2,2>*V(3,3)-V<  1,  3  )*V  (  2,  2  )»V<  3,  1  )+V(  1,2)*  V(2,3)»V(3 

00250C 

' 

Ii1)-VI1,2!»V(2,1)»V(3,3MV(1,?)»V(3,  2)»V(2,1)-V(1,1)»V(3,2)»V(2,3) 

00251C 

GO  FC  61 

002520 

21 

PDEL=DEL(L) 

0C2550 

GO  TC  25 

0C25UC 

22 

IF(AA)23,23,2U 

00255C 

23 

PDEL  =  DEL(L-1  ) 

0C2560 

GO  TO  25 

002570 

2k 

PDEL=CELC 

002580 

25 

IF(PCEL*DEL(L) 126,26,27 

002590 

26 

IF(AA)28,28,57 

002600 

21 

BB  =  0 

002610 

IF(AA)3b,35,52 

OC262C 

28 

AA=1.0 

002630 

B6=l  .0 

0C26U0 

DELC=LEL(L-1  ) 

002650 

CfEGAC  =  C^EGA(L-1  ) 

00266C 

GO  TC  59 

002670 

2900MEGAMtM)=OK'E 

0C268C 

1 

IABSF(CEL(L)  )  ) 

002690 

GC  TC  32 

0027CC 

30 

CMEGAIL+1  )=CMEGA(I.  )-.9«UCf 

00271C 

D0M=.5»LCI*  OC27JC 

31  L  =  L+1  .  0027-.0 
IF( JK)  100C,  10CO,  1002  0C27U0 

32  BM=M  OC27  50 
IT  (HN-eV)3000,3CCC,  ii  00276C 

33  M  =  M+1  002770 
IF (N-313UC, it ,3U  0C27PC 

34  AD0M=(0MEGAM(C-1  )-0MEGAM(M-2  ) )  •.  5:  002790 
CR=ACrw/bC.  002CCC 

340  AA=0.0  OC281C 

BP*0  0C28':0 

COC  0C2G3C 

CG  =  C.O  0022UC 

DOK=ADCM  0C23'jC 

DELC.  =  0  0C2S60 

GMEGML«-l)  =  CMEGAy(M-l)+2.»Cfl  0  026  7C 

CHEGAC  =  0MEGA(L*1  I  002JirC 

GO  TC  31  00289C 

35  IFUNECA(L)-CyEGAG)  55,60,  6C  0C2900 

50  0MEGA(L+1 )=CMEGA(L)-(nM:GA(L l-CMfCAC )/?.C  002910 
IF(CCfcGA(L)-OfEGA(L  +  l )-GR)51  ,51, 31  0  02  92C 

51  COI.O  002930 
GO  TC  31  002940 

52  OELC=UEL(L)  OC29bO 
CKEGAC=OMEGA( L)  ,  00296C 

53  0MEGA(L+1 1 =CMEGA ( L ) +DOV  002970 
GO  10  51  0C298C 

54  IF(LELC»DEL(L)).?9,29,55  0C299C 
550CMEGAMMI=0HEGA(L»  +  (CMEGA(L-1 ) -GKEGA ( I ) ) *ABSF ( DEL ( L ) ) / { A BSF ( DEL < L )   0C3C00 

1 )+ABSF(CEL(L-1  )  )  )  0C3C10 

GO  TC  32  OC3C20 

56  IF(CNEGA(L)-OKEGAC)21,21  ,22  0C303C 

57  IF(BB)5e,58,5G  003C40 

58  PR=1.0  OC3C50 

59  IF(0VEGA(L)-0I-EGA(L-1  )-CR)29,29,  ^0  0C3060 

60  M=M-1  003C7C 
GO  TC  3G00  0C3C80 

61  IF(GG)70,70,72  OC3090 

70  IF(DEHL) 177,71,77  003100 

71  GG=1.0  003110 
OMEGA (L+l J=OMEGA(L)-CR  OC3120 
GO  TC  31  OC3130 

72  IF(GG-1. 0)73,73,7U  0031U0 

73  GG=2.0  003150 
HOEL=CEL(L)  003160 
OMEGA (L+l  )=CMEGA(L)+2.»CR  0C3170 

.  GO  TC  31  OC3180 

74  TOEL=DEL(L>»HCEL  003190 
IF(TCEL)75, 76,76  0032CO 

75  0MEGAM(M)=CMEGA(L-2)  003210 
GO  TC  33  0C3220 

76  KK=1  003230 
>=M-1  0C3240 
GC  IC  30CC  003250 


77  IF(CC)56,56, 54  0C326C 

C  0C3270 

C      OUTPUT  0C328C 

C  0032^0 

3000  IF(NM)4C0,40O,4Ci  0C3iCC 

400  PRIM  UC1  0C3310 

U01  FORMAT  (  1H1//  )  003320 

PRIM  4C2  00333C 

4020FCRMAT  (15H  PROGRAM'  VI  P  IPC,  35Xt  17HG.  E.  FINK   NHA- 1  ,  55X,  1  3HFfcBRLAR   0C3340 

1Y  1964//,  10X,  TOSH  MODAL  FREQUENCIES  OF  A  PLANAR  PIPING   0C3350 

<.     SYMEM  ARE  DETERMINED  pv  AN  ITERATIVE  PRCCTCURE  USING  T,HE  /29H  ME   OC3360 

3THCL  CF  TRANSFER  M  A TP I C E S . / / /  1  1  9h  •••••••«•••••••«•*•«•••»•••»•«•»   00337C 

^••••••«b««*«»»«««««««<««««»*»»«««««« ••*«>••*•*•<«•**««•« ««••••■>«*   003380 

i //)  0C3390 

403  IF(NM»405F,  4059,  4C4  0C34CC 

404  PRINT  4C5  003410 

405  FORMAT  (  1  HI / /  )  0C3420 
4059  NM=i  0034J0 
4059  IF(  JK)4(.6C,4060,IjC65  0C344C 
4C60  PRINT  4070  JJ  003460 
40700FCRMAT  (  5CX  ,  3HPR0PI  EH     I  2  ,  /  /  t  5  X  t  M  7H  IM  PLANE  MODE  FREQUENCIES  (RAIMA   003460 

INS  PER  SECOND)    //45X,  Ut-VOCfc,  15X,  10H  FREQUENCY  /)  003*470 

GO  IC  4079  00349C 

4065  PRINT  1+075  JJ  003490 

4C750FORMAT  (SCX,  3HPR<^LEM  I  2  ,  /  /  3  :.  X  ,  5  1  HOOT  OF  PCANh  MCCF  FREQUENCIES  (RA   0C3500 

1DIANS  PER  StCCNC)       /  /45X  ,  UHMODEi  1  5X  ,  10H  FREQUENCY  /)        003510 

4079  MH  =  M  00352C 

PRINT  408  tM,CML"GAK(M),M=l,MH)  003530 

408  FORMAT  ( 46X, I  2 , 1 7X , F 1 4 . G  /)  003540 
IF(KK)4C80, 4090,4080  0C355C 

4CP0  PRINT  4035  0O356C 

40e50FORMAT  (  1  OX, 69FPRCHL EM  TERMINATED  CUE  TO  SIGNIFICANT  FIGURE  LIMITAT   00357C 

HON  CF  COMPUTER.      /)  003580 

4090  IF( AMYK)412,409,4l2  003590 

409  PRINT  U 1 1  003600 
4110F0RKAT  (16X.53H               A  LUMPED  MASS  APPROACH  WAS  EMPLOYED,  A   003610 

1N0  )  0C3620 

GO  TO  414  003630 

412  PRINT  413  003640 


412  PRINT  413  003640 

413  FORMAT  (24X.45HA  DISTRIBUTEE  MASS  APPROACH  WAS  EMPLOYED  AND       )    003650 

414  IF(SD)415,415,420  003660 

415  IF ( R I  1416,416,418 

416  PRINT  417 
4170F0RMAT  (2X,fc7hTHE  EFFECTS  CF  SHEAR  JEFLECTION  AND  ROTARY  INERTIA  W 

1ERE  CONSIDERED.  /) 
GO  TC  425 
418  PRINT  M9 

4190F0RMAT  (3X.63HTHE  EFFECTS  CF  SHEAR  DEFLECTION  WERE  CONSIDERED  WHIL 
IE  THOSE  OF      /3X, 35HROTATICNAL  INERTIA  WERE  NEGLECTEC.        /) 
GO  TC  425 

420  IF(RI)421 ,421  ,423 

421  PRINT  422 
4220F0RMAT  (3X,7EFTFE  EFFECTS  CF  ROTATIONAL  INERTIA  WERE  CONSIDERED  WH 

1ILE  THCSF  OF  SHEAR  / 3 X , 2 7HD E FLEC T I  ON  WERE  NEGLECTED.    /I 


C-C  IC  1*25 
423  PRIM  U2h 
l*2«K)FCRyAl  (^A.treTHF    EFFECTS 

IRC    NEGLECTED.     /) 
U25    PRINT    U^C 
l*i+0CFCRPAT  (  11°H    «*.«««.«*»*« 


ROIAr-!Y     INEMia 


IFICDABC-1.  )U51  I  ,  "4  5  2C,  «*i4C  1 
UUOl     PRINT    UiiC2 

l»l»020FnRFAT(;>5H    SECTION    LENGTH    ANL    SUPSE: T  ION    DAT  A/2  3X,  1UHSE:  TION    NLKRE 
1R,6X,;5HLENGTF    OF     S  tC  T  ION  (  I  NCHES  )  ,  5<  ,  2  1HNUMFJ  Eu    OF     SECTIONS  ) 

PRIM    km     {  I  ,TL(  I  )  ,Z(  I  ),  1=  1  ,NS,   ) 
UU3    FORMAT    [29X,  I  3, 23X, F7 . 2 , 2 1 X , FU.O ) 
LC  =  L 

PRIM     UM 
l*l*10FCRKAT     (11H    GRAPH    DATA    /  2  '  X  ,  1 6  H  I  T  LRA  T  I  CN    NUM  P  ER  ,  K  X  ,  5  CHFR  CCUcNC  Y     (R 
1ADIANS    PEP     SECOND)  ,  3X,*r>F  VALUC    Cc     FRlGUCN'CY    DETERMINANT  ) 

PRIM    UUU     (L.CPEGACLl ,CEL(L),L=1,LM) 
UUU    FCRf'AT     (29X,I3,23X,  F11  .6, 18X.E15.8) 
'♦520    PRINT    UU2 
1*1*2    FHR.VAI  (  IX,  10HINPLT     [ATA     /) 
PRINT    U5  21 
4521  OF  OR* AT  (  1 X , 9HLCVPCNENT , 4 X , RFl  I  AVE  TFR,  5X,  14HWALL     THICKNESS     ,5X,6FL£N 
IGTH,5X,9HRACIUS    OF, 9X , SH INCLUDED , 5X,  7HCLNS I T Y , 5X , 7HE LAST IC , 5X 1 5HSH 
2EAR    /c.7x,9HCUPVATLRl,SX,  12FANGLE    CF    ARC  ,  1  6X  ,  7HMC  DLLU  S  ,  5X  ,  7HI"CDULUS 
i    /) 
Dr    U525    1=1, NS 
IF  (PI  I  )  )  4 52 4, 4  5  2  4,452  3 
4524    PRINT    4522     I  ,C(  I  )  ,LT{  I  )  ,TL(  I  )  ,RHC'(  I)  ,  TFMR(  I  >,BAKU(  I  )  ,E(  I  )  ,G<  I  ) 
4  522  0F0RKAT     (  3X  j  I3,8X,F7  .3,  10X,  F6  .3,  7X,  F3  .  3,  4X,  F8.3,  9X,F6.2,P  X,F7.2,**X, 
1EF.2.4X, te.2 ) 
GO    TC    4525 
4523    SH=1. 
U525    CCNTINUL 

IF(Sh)U53G,4S30,4526 
45?6    PRINT    U532 
4532    FORMAT     (// 1 X , 7HH ANGERS    /) 

PRINT    U527 
l*5270FCRKAT     (  2X,  9HCCf  PCNENT  ,  1  4  X  ,  ;HCLX  ,  1  4X  ,  3HCLY,  1  4X  ,  3HC  LZ  ,  14X  ,  3HCTX,  1  U  X 

1 ,3HCTY, 14X,5HCTZ       /) 
1*528    CC    U530    1  =  1  ,NS 

IF(P(  I  )  )  1+530, U530,«*  529 

4529  PRINT    1*531     I  ,  CL X (  I )  ,C LY (  I  ) , CLZ ( I  ) ,CT X <  I ) ,CT Y  (  I  )  ,C TZ (  I  ) 

1*531     FCRMAT(iX,I3,15X,F1C.:,7X,FlC.2,7X,FlC.2,7X,F10.2,7X,F10.2,7X,FlC. 
12     ) 

4530  CONTINUE 
PRIM  Ul*8 

1*1*8  FORMAT  (7/20H  BOUNDARY  CCNCITICNS  ) 

PRINT  449  (  SVK  I  ,  1  )  1  =  1  ,6  ) 
1*49  FCRKAT  (5H  SVI=6F3.0) 

PRINT  U51  (SVE( 1,1)  1=1, 6) 
1*51  FORMAT  (5H  SVE=oF3.0) 
IF(NBR)451  1  ,1*51  1  ,  U5C2 
U502  IF(JK)U505, 4505, 4515 
4505  PRINT  «*510(  (SVIP(  I,  Jl  1=1,6)  J=2,NNN) 


0C38  <0 
0C3'JU0 
O03c?5C 
00  3860 
0C3870 
003P.«C 
0  0  3  8  9C 
OC39C0 
0C3910 
0C391C 
0039:SC 
00394C 
0C395C 
0059tC 
0  0  3  9  /  C 
00398C 
0C3990 
0C14C0C 
00401C 

ocucro 

004030 
004040 
001*050 
0C4G60 
00UC70 
001*080 
00UC90 
0C4100 
0041 10 
00U120 
004  1 30 

004140 

00U150 
00U16C 
0CU170 
0C4180 
0C4190 
004200 
001*210 
0C42^C 
00U230 
00U235 
0C4240 
0C4250 
001*260 
00U270 
00U2SC 
004290 
00U3C0 
00431C 
004320 
004330 


1*510  FCRVM  UF3.0)  00U5U0 

GC  TC  U511  OCU35C 

1*515  PRINT  i4b  1  7  (  (  bVCP  (  1 1  J  >  1=1,6)  J  =  2,NM)  004560 

4517  FPR'fAT  (<^F3.G)  004370 

U511  CONIIMJ:  OCUiBC 

nr  kkb   iw*i,i*h  0CU390 

i*U6  riMECAMN)=0.0  0044CC 

nr  447  l=i,lm  ocuuio 

CPtGA(U=C0  004420 

447    DEL(U=0.0  0044<0 

DELC=C.C  004440 

CfEGAC'0.0  00445C 

IF(MN)453,453,452  004460 

1*52    MN«0  004470 

L=i  oc44ec 

jwt=  i  00449C 

AA*0.0  0C45C0 

BB=G.C  00U510 

CC  =  0.0  00«<520 

GG=0.0  00Ub30 

KK»C  0C454C 

ADOP*hOOM  0CU550 

CR  =  f>R  .  004560 

0MEGA(1>=0MLGM  004:>70 

OMEGAC»OMEGAl    '  0045P0 

DOM=ADOM  004590 

GO  TC  908  004600 

U53  JJ»JJ+1  004610 

DO  454  1=1, NS  004620 

P(I)=0.  004630 

454  Z(I )«C.C  004640 

GO  TO  U50  0C4650 

tNO  00U660 

C  004670 

C  004600 

SUBROUTINE  SUeSEC  I  RHP, TFETA , C, TL,  FN  ,  Z )  004690 

RH0U=ABSF(RH0)  004700 

IF(HV-5.0>8,0,7  0014710 

7  HH»S.O  001*720 

8  RA»TL/D  004750 
IFIRHOD-0. )3, 2,4  004740 

3    IFIRA-1. >5,5,6  004750 

5  Z»0.  004760 
RETIRN  0C4770 

6  IFIRA-3.)  15,15,16  004780 

15  Z«2.0  001*790 
•  RETURN  00U8C0 

16  IFIRA-6. 117,17,19  004810 

17  Z»3.0  004820 
RFTURN  004830 

18  IFIHM-3.11, 1,2  004840 

1  Z»6.C  ooueso 
RETURN  004860 

2  Z=2.0»HK  004670 

D-10 


RETURN  0048"0 

k    RO=RHCD/n  00.U89C 

IFCRL-1.  )19, 19,23  0049CC 

19  1=0  0C491C 

RtTLRN  OC4920 

23  IF(.RL-3.  )?4,24,25.  0CU93C 

24  [FITHETA-. 1*4)26, 26, 27  004940 

26  2=0.  0C4950 
RETURN  0C4960 

27  IF{  THITA-1  .37)26,  13,29  0C49fC 

28  Z=2.0  00491C 
RETURN  O0U99O 

29  Z  =  3.C  0C5CC0 
RETURN  00SC1C 

25  IF(RC-6.)30,3C,i1  0C5020 

30  IF(ThETA-.44)22,32.,33  0C5040 

32  Z  =  2.C  0C504C 
RETURN  00505C 

33  IFUI-LTA-1.C,7  ) '4,34,35  005C60 

34  Z=3.0  005070 
RETLRN  005COC 

35  Z=4.C  005C9C 
RETURN  0051C0 

31  1F( THETA-.44) 36, 36, i  005110 

36  Z=2.0  '  00512C 
RETURN  005130 
ENO  005U0 

00515C 
0C5160 

SUBROUTINE  DISTN  ( AMU , A  I Y ,OPEGA, 1 L ,3 , C, 0  I  ,C, SD, R 1 , E , FT AC ,U  )  00517C 

DIMENSION  U(6,6)  005100 

ASA=<3.1415  926  54/U.0)«(C**2-CI»«2)  005190 

AREA=ARA*FFAC  OC52CO 

PE=TL»GMCGA*SCRTF(  AMU/(E«AR/s  )  )  0052  10 

PU=R«AMU»TL«»U«CKEGA»*2  005220 

IF(RI)1,1,2  005230 

1  T=R*ANU«(  AIY»C;GA*TL)»»2  0C5240 
GO  TC  3  005250 

2  T=0.0  0C5260 

3  IF(SC)4,4,5  005270 

4  S=(AKU»(CfEGA»TL)««2)/(G«AR;A)  005290 
GO  TC  6  005290 

5  S=0.0  005300 

6  SPT=S+T  005^10 
VV  =  SCRTF(PU+( S-I  )«*2/u.0 )  00532C 
AL1=SCRTF( VV-.5«SPT )  005330 
AL2=SCRTF< VV+.5»SPT)  005340 
CL=1  .0/1 AL1*«2  +  AL2»»2 )  005350 
CHL1=CCSH(AL1 )  00536C 
SHLl=SINH(ALl )  005370 
CL2=CCSF( AL2)  0C5380 
SL2=SINF(AL2)  005390 
SPE=SINF(RE)  0054C0 
CBE  =  CCSF(f<E)  005410 


CCO=CL«UL2»«2«OL1+ALU»2*CL2)  OOSUTC 

CC1=CL*UAL«1««2»5I-U)/AL1+(AL1**:*SL2)/4L2)  0  054/50 

CC2=CL»(CHI_1-CL2)  0C5440 

C03=CLM SHL.1/AL.1-SL2/AL2)  00  54  5C 

DC  7.  J*l,6  0  0  546C 

l.'C  T    K=l  ,t  003U70 

UUtK)=O.C  0C54T0 

L(1,1>=LHE  0C549C 

U(l,6)  =  TL«SPE/(Ct»A-RA*E)  OC5r.CC 

U(2,2)=CCC-S»CC.'  OC5S1C 

U  (2,  3  )  =T  L*  (  f.C  1 -SPT  «CO  7,  )  0C5S20 

U(2,U)=R»CQ2«TL««2                                      J  00553G 

U(2,5}=(R«TL»*3/P4).((P4+S««:)»CL3-S«CC1)  0  05  S4C 

U(  3,2)=PU«CC3/TL  0C555C 

1.(3,  ;)=CCC-T«CCZ  005d60 

LM3,4)=TL*R«(CG1-T»C0?J  005570 

U(3,5)=U(2,4)  0055>n0 

P o  =  1  .  C / 1'  0C5590 

U(4,2)=P4eRR»CC2/TL«v  005600 

U(4,3)=RR»((P4+r««2)«CC3-T»CCl)/TL  005610 

U(Uf4)=L>(  3,  3)  OC5620 

U(4,5)=l(2,3)  005630 

U(5,2)=P4»RR« (CC1-S*C03>/TL««*  005640 

L(5,3)=L(4,2)  0C565C 

U(5,4)=U(3,2).    '  00566C 

U(5,5)=b(2,I )  0C5670 

U(6,l)=-AKJ«TL«CN'EGA»*2«SFE/HF  005690 

U(6,6)=CBE  00569C 

RETURN  0057CC 

END  005710 

00572C 
00573C 
SUBROUTINE  CISTVC  ( AND, A  I  X , A 1 Y , P^EGA , TL , AJT , R , 0 , D I , G, SO,  R I , FF AC , L )   0C5740 

DIMENSION  U(6  ,6)  005750 

k*l./(G*"AJT)  005760 

AREA=(3.1U159265U/U.0)*(r.»«2-DI*«2I*FFAC  005770 

BF«*SCR1F(AVU  •  (  T  L   » A  I  X »CMEG A ) •» 2« k )  005780 

P4=R«AMU    *  TL  ««4*0yEG/S»«-  005790 

IF(RI)2,2,3  0058C0 

r*R*AMU*tAIY»CPEGA»TL)«*2  005810 

GO  TC  4  005820 

T=0.0  005830 

IF(SG)5,5,6  005840 

S=(ANU»(CVEGA»TL)*»2)/(G«ARLA)  005S50 

GC  TC  7  OC5060 

S=0.0  0C5870 

VV=*SCRTF(P4+  (S-T)'*«2/U.C)  005890 

SPT=S+r  005R90 

AL1=*SIPTF(VV-.5»SPT)  0059C0 

AL2=+SGRTF( VV*. j»SPT )  005910 

CL=  1  .0/( AL1»»2  +  AL2»«2 )  00592C 

CHL1=CCSH(AL1 I  005930 

SHL1=SINH( AL1 )  O059U0 

Cl2=C0SF( AL2)  0C5950 


SL2  =  SIW  UL2  ) 
SBE=SlNF(tEI 

CPE-CCSr(rE ) 
.       CC0  =  CL«(AL2«»:*CHI   1     ♦ 

C01=CL*(  (  AL2»*2»Sf-Ll)i 

C02=CL»(CHL1-CL?) 

C03=CL«(SHL1/AL1-SL2/All) 

DH  1  J=l,d 

DO  1  K=l,6 
1  U(J,K)=C.C 

uti , 1)=cpe 

U(,l  ,2>=TL»W«SRE/BE 

U('2iU=-AHJ   »  TL   •AIX*«:«SPE/Rl 

U(2,2)=CBfc 

IK3,3)=CCC-S»C02 

U(3,U>=  TL  *(CC1-St*r*COi) 

U( 3t5)=R«CC2«TL««2 

U(3,6)  =  (R«rL»«?/Pit)«(-S*CCl-MPi4+:.«»2 

U{Ut3)=PU«CC_">/TL 

U(u,4)  =CCC-T*CCL 

U(U,5)=  TL»R»(CC 1-T»C03) 

U{U,6)*U{3,5) 

RR=1.C/R 

U(5*3)=P4*Rr<«CC2/TL»«Z 

U(5,t»)=RRM-T«CCl*(P»»«r**Z)«CO*)/TL 

U I S ,  5 )  =U  ( U  f  U  ) 

U(5,6)=U(3,U) 

U(6,,3)=PU»RR.(Cri-S«cn3)/TL«»3 

U(6,li*  l=U(c,  3) 

U(6,b)=L(4f :) 

U16,6)*U(3, «) 

PETLRN 

END 


(All «»2*st:  )/AL.  ) 


SUBROUTINE  SFIEIC  (fl.TL, 

OIKtNSION  A(6,6) 

ARA=(3.1U159?65U/a.C)«(r.< 

AREA=ARA»FFAC 

UL  =  11/1 

DO  1  J  =  2,  5 
1  A(1,J)=0.0 

A (1  ,  1  )  =  1 . 0 

AM  ,  6)=DL/(  l«ARA) 

A (2,  1  )=G.0 

A  (  2  ,  2  >  =  1  .  0 
.  At2,5)=DL 

A(2,U)=DL« 


P  ,  /  ,  G  ,  S I  ,D,  E.FFAf.,  A) 


IF(SU)2,2,3 

A(2,c;  )=l  !_•• 

.»R/fc.C- 

UL/(G*ARLAJ 

GC     IC    U 

A(2,5»=UL»» 

«R/6.r 

A(2,6)=C.O 

A(3,  1  )=C.O 

0C5960 
005970 
00598C 
0P599C 

nc60Co 

006010 
0C6LVC 
00603C 
006040 
0C6Cd0 
0C6C6C 
0C6G70 
006080 
00609C 
006100 
0061 10 
0C612C 
00613C 
0061UC 
006150 
0C6160 
006170 
00618C 
00619C 
006200 
006210 
00622C 
006230 
006240 
0062^0 
006260 
006^70 
O062B0 
00620C 
006300 
006310 
006320 
0  G  6  3  3  0 
006340 
006350 
006360 
00637C 
0063P0 
006390 
006400 
00641C 
006420 
006430 
006440 
0C645C 
006460 
0C647C 
00648C 
00649C 


A (3,2)  =0.0 
A  (  3  ,  5  )  -  1  .  0 
A(3,i»)=uL«R 
A  ( 3 ,  5  J  =  A  ( Z , 
A(3 ,61=0.0 
DO    5    J  =1,6 

a(u, j)=c.c 

A (5, J)=L.C 
A (6, J) =0.0 

A  (  U  ,  h  )  =  1  .  0 
A(U,5)=I>L 
A(5,S)=1  .C 
A  (  6  ,  o  )  =  1  .  0 

RCTl;R\ 
END 


SURKCU 
DIMLNS 
AREA=( 
UL=TL/ 
DC 

r,c 

A( J,K) 
A{1,  1) 
A  {  2  ,  2  ) 
A(3,3> 
A(U 


M5, 


riNE  SFIELf  {  CI  ,  n  ,  A  J  I  ,R,Z,G,SlJPl,l),FFAC,  A) 

ION  A (6, 6) 

3.1U15T26  5i4/U.O)»<r.«*;-ri**2)«FFAC 

z 

=  1,6 
=  1,6 
=  0.0 
=  1  .0 
=  1.0 
=  1  .0 
-1.0 
=  1  .0 
=  1  .0 

. )7, 7, R 
=0L»«i«R/6.C-CL/(G«ARLA) 


=  IA**2«R/2.C 
=UL»R 
M3,5  ) 


A{6,6) 
IFISCR 
A(3,6) 

cc  rc 

A(3,6) 
AM, 2) 
A(3,M 
A(3,5) 
A(U,5) 
A(U,6) 
A(5,6) 
RETURN 
END 


SUBROUTINE  CFIELU  (R, 

DIMENSION  A (6, 6) 

PHI=THETA/Z 

RD=AeSF( RHO) 

IF(KHf ) 1  ,  1,2 

V  =  -l  .0 

GO  TO  3 

V=l  .0 

CP=LCSF(PFI ) 

SP=SINF(PHI ) 


Z,G,SD,RHC, THETA,n,OI,E,A) 


oo6l:oo 

0  0  6  5  1  C 
0065 iO 

CC6-j:C 
0C6^UC 
0  0  6  5  5  C 
0C6560 
OC6  5  70 
0C6bP0 
00659C 
0C66C0 
0C6610 

oc6t:c 

00664C 

006640 
006o jO 
00666C 
006670 
0C66°C 
0C6690 
0067C0 
006710 
0C6720 
006730 
0C67UC 
006750 
0C676O 
0C6770 
0C678C 
006790 
OC66C0 
006S1C 
006E2C 
006830 
006CU0 
0C6850 
006C6C 
006870 
0C6830 
006890 
0069C0 
006910 
0C6920 
0C69  3C 
006940 
006950 
006960 
006970 
0069«C 
0C6990 
007000 
00701C 
007C20 
0C7C30 


IT  U  J=l  ,6 

cr.  «  k  =  i  ,6 

A ( J,K)=(  .( 

&  <  1  ,  1  )=CP 

A( 1  ,2J=-V*SP 

A( 1 ,3)=-V»RC« ( 1 .O-CP) 

A{  1  ,U)=-V*RC««2«R.»(PHI-SPI 

AR=3.141592£5U*(C**2-EI*»2)M.C 

U=1.0/(t#AR) 

F3=(SP-PHI#CP)*.5 

F1=IPHI»CP+SP)«;5 

F5=(2.0-2.0»CP-PHI»SP)«.= 

F6= ( 2 . 0*PHI +PHI »CP-3 . 0«SP ) * .  5 

A  (  1,  5)=V»(R  C  «U«PH  I  «SP  «.:  -P»l  J»T 

An,6)=RD«W*Fl+R*F6*RD«*3 

A ( 2  » 1 l=V»SP 

A(2,2)=CP 

A(2,3)=I'D«SP 

A(2tU )=R«RD*»2« ( 1 .O-CP  ) 

A(2,5)=RP*F  5*U*R«RC*»2) 

A12,6)=A(1,5) 

A  (  3  ,  3  )  =  1  .  0 

A(3,U )=RC»R*PhI 

A(3,5  )=A(2,H) 

A(3,6)=A(1,i4) 

A  { M  i *» I  =  1 . 0 

A(U,5)=M2,7) 

A(U,6>=A< ]fi) 

A(5,5i)=LP 

A(5,dl=A(  1,2) 

A{6,5)=A(2,T) 

A(6,6)=CP 

RCTLRN 

END 


SUBROUTINE  CFIELC  ( AJT , P , I , C , SCR  I , RHC , TH ETA, D,P I , FF AC, A) 

DIMtNSICN   AI6.6) 

PHI=THETA/Z 

RHGI,  =  RHC 

IF(RHO-O.O) 11,11,12 

11  V  =  -l  .C 
GC    1C    1 

12  V=1.C 

l   CP  =  ccsF<r»h'i  ) 

SP=SINF(PHI ) 
W=l.C/(     C    »AJ1 ) 
F1=(PHI»CP+SP)/2.C 
F3=(SP-CP«PHI  )/?.(' 
F5=(2.0-2.0«CP-PHI*SP)/2.C 
F6*(2.O*PHI  +  PH*CP-3.O*SP)/2.0 
RHC=ARSr (RHC) 
,    DO    „     J=l  ,6 


CO 


K=l, 


007CUC 
0C7C5C 
CCTCrO 
007C70 
O07C- C 
0C7C9C 
0071C0 

o  c  7 ;  i  o 
00712C 

OC7  1 ;0 
C071U0 
OC71oC 
0  C  7  1  r-  C 
00717C 

co7; . o 

00719C 

0C7^C0 
0  C  7  "  1  0 

oc7:,:c 

OC7_rsO 
0C7240 
OC7250 
007260 
007^70 
OC723C 
0C7290 
007^00 
007/S1C 
0C7i~C 
0073.^0 
007<U0 
007i5C 
0C7360 
007^70 
0075S0 
0C739C 
0C7UG0 
007U10 
0C7U2C 
007UiO 
OC7UH0 
007U50 
0C746C 
007470 
0C7420 
0C7U90 
007i:00 
007d)0 
0C75P0 
0C7b!C 
0075UO 
007550 
OC756C 
007-70 


2  A(J,K)=C.C  0C7S:iC 

IFIbDRI ) 15, 13.1M  007590 

13  AP.E.*=(3.1Ulb9265U/t».0)»(C**:—  DI*«2)»FFAC  007tCC 

Sn»RHC»THETA/(G*Z»AREA)  007fclC 

GC  TC  lb  OC762C 

1U  SD*0.C            •  C07tiC 

15  A(  1  ,  1  )=LP  007CUC 

A(l,2)=(k»RHG«Fl)-(RHO«F3*R)  DC76',0 

A(1,U)=V«SP  0C7660 

A<1,5)=V«(W+R)»(RhO«PHI»SP)/2.0  0C7t7C 

A(l,fc)=V«(W+R)«(RHC»«2)«F2                              -  007620 

A(2,2)=tP  00769C 

A(2,5>=V«SP  0077oC 

A(2,6)=V»»HC»( 1 .G-CP)  0C771C 

A(3,l)=-A(2,6)  00772C 

A(3,2)=-1.0»A( 1,6)  OC775C 

A(3o)  =  l.C  O077HC 

A(3,U)=RHC*SP  0C77j0 

A(3,5)=((R»<RH0««2)*PHI»SP>/2.C)-W*RhC»*2»F?  007760 

A(  3,6)  =  (R»RHO«*i»F3)-U»F6»RHC»»  -  ) -S 0  OC7  7  70 

A(U, 1 )=-V»SP  007780 

A(U,2)=-V»(W*R)«RF0»PHI»SP/.:.0  007790 

A(U,t)=CP  007800 

A(U,5)=(R«RHG»F1  )-(W»RFC«F3 )  0C7810 

A(i4  ,6)=A(3,5)  0C7820 

A(5,2)=-V»SP  0C7B30 

A(5,5)=CP  0078UC 

A(5,6)=A(3,U)  OC7850 

At6,6)=l.C  007860 

RHO*RHCC                        .  007670 

RETURN  0078°.0 

END  007890 

007900 
007910 

SUeRCUTINE  POINT  ( A^U, A  I Y, OP  EGA, FL ,Z , R I , B )  007920 

DIMENSION  B(6,6)  0079.30 

AM=AMU«TL/( Z-l.  )  0C79U0 

DC  1  J=l,6  0079^0 

DC  1  K=l,6  0C7960 

1  B(J,K)=0.C  0C7970 
B(l,l)=1.0  OC79P0 
P(2,2)=1.C  007990 
H(3,5)=1.C  008CCO 
IF{RI)2,2,3  0C8C10 

2  R(U,3 )=-Ay*(AIY«CVEGA)»«2  008020 

3  B(U,U)=1.0  008030 
B(5,2)=AK»OVEGA««2  0080*«0 
P(5,5)=1.0  008050 
P(6, 1 )=-P(5,2 )  00806C 
B(6,6)=1.0  008070 
RETURN  OC808C 
LND  008090 

OC81C0 
008110 


SUBRCUTINF  POINC  (  A  VU  ,  A  I  X  ,  A  I  Y  ,  ON  '  G  A,  T  L  ,  I  ,  SD*  I  ,  P  )  OC8120 

DIMENSION  P. (fr,6)  C081 5C 

AN  =  AMI«TL/(Z-1.C  )  OCB1UC 

DC  1  J=l ,6  OC31 dO 

DO  1  K=1,6  OC8UO 

1  P  (  J,K  )=C.C  00  9  i  ,'C 

P(1,1)=1.C  0081^0 

B(2,l)=-Ay*(AlX»nwECA)»»2  0C8190 

P(2,2)=1.0         •  oc8:co 

P(3,3)=1.C  0C8210 

B(U,«4)  =  1.0  0C8220 

IF(SDRI-0. )9,9, lc'  0C8230 

9  P(5,4)=-AM« ( AIY«CNEGA) »«2  008^40 

GO  TC  11  0C8230 

10  6(5, M=0  0C81O0 

11  E(5,5)=1.C  008270 
P(6,3)=AM«CN!"CA««:  008.^0 
P(6, 61=1.0  00829C 
RETURN  0C8300 
ENC  0  08.1- 1C 

00832C 
008330 

SUBROUTINE  NATMUL  (A,R,Z,U)  0C83U0 

DIMENSION  A(6  ,6)  ,F(6, 6),C (6,6),D(6,5  > ,U(6,6)  0C8350 

DC  1  J=l,6        '  00836C 

DC  1  K=l,d  008370 

C(J,K>=0.C  0C8380 

DO  1  L=l,6  0C8J90 

1  C(  J,K  )=(.(  J,K)+B(  J,L  )»A(L,K)  008400 
N'  =  Z-1.8  008410 
DC  Z  J=l,6  0G8420 
DO  2  K=l,6  008*430 

2  U(J,K)=C(J,K)  008440 
IF( Z-2.0)6,6, 3  008450 

3  DO  5  M=l ,N  008460 
DO  4  J=l  ,6  0C8470 
CO  U  K=l ,6  OO840C 
D(JfK)=0.0  0084<90 
DO  4  L=l,6  008500 

008510 
008520 

DC  b  K=l,6  008530 

008540 
00855C 

DC  7  K=1,6  0C8S60 

D(J,K)=O.C  008570 

'  DO  7  L=l  ,6  008580 

7  D(  J,K1=D(  J,K)*A(  J,L)»U(L  ,K)  0085*30 
DC  j  J=l,6  008600 
DO  B  K=l,6  0C8610 

8  U(J,K)=D(J,K)  008620 
RETURN  00863C 
END  008640 

008650 


(  AMJ,  TL,CM:GA,A  [Y,l  I  ,L  ) 


SUBROUTINE  RIGI 

DIMENSION  U(6,6) 
AV=AM)«TL 

nc  i  j =  1,6 

DC  )  K=1,6 

L( J,K)=C.C 

U  (  1  ,  1  )  =  1  .  0 

U(6 , 1 |a-AV»CMEGA»»2- 

U  (  6  ,  6  )  =  1  .  C 

U  ( 2  ,  Z  )  =  1  .  0 

U(2,  31-TL 

U  (  3 ,  3  )  =  1  .  C 

U(«t»2)=AMTL»CMEGA»«2/2.C 

IF(*I)2,2,3 

U(U,.5)=AM«Ci*ECA«»2»(TL«»2/6.0 

gc  rc  »* 


U(U  ,  »  )=AP«{  CM  EGA 
U  ( 4 1 4  )  *  1  .  0 
U  (  U  ,  5  )  =  I  L 
U(5,2)=-U(6,l ) 
U  (  5  ,  3  )  =1 ( 14 ,  2  ) 
U(5t5)-1.0 
RETURN 
END 


TL  )  » 


2/6. C 


til,  1  ) 

U(2 

U(2 

U(3,3) 

U(3,U) 

U(U,U) 

U(5,3> 

I  F  (  SIJR 

U(5,U) 

GC  IC 

U(5,iO 

U(5,5) 

U  ( 5  , 6  ) 

U(6,3) 

U(6,U) 

LM6.6) 

RETGRN 

enii 


UTlNt  RIGIC  (  AMU,  T L  ,  A  I  X , CKtGA ,A  IY,  AJT,  SCPI  ,U) 
SIGN'  U(6,6) 
J=l,6 
K=l  ,6 
1=0.0 
U*TL 
=  1  .0 

MU»TL«( aix«lwe:a)*«2 

=  1  .0 
=  1.0 
=  TL 

1.0 
=AM«TL»CMEGA»«?/? .C 
I-C.)  3C,70,31 

AM«OMECA»»?«(Tl«*2/6.0-AIY«*2  ) 


=  AM«C:'EGA 
U(D,3) 
1.0 


.0) 


SUBr<CUTINE  STIFCC  ( THET A  ,RKC ,U ) 


ocet6  0 

0C8670 
008680 
00369C 
0087CC 

ocsnc 

008720 
0  0  8  7  3  C 
OC87.C 
OC8750 
0C87f-.C 
0087^0 
008730 
0C8790 
008600 
008° iC 
008820 
0088  '0 

ocaaiio 

ooe  n 5 c 

0C886C 
008670 
008680 
00889C 
0C89C0 
0C8910 
0C892C 
0C89-50 
008940 
0089'jO 
00896C 
0C897C 
0C89oC 
008990 
00900C 
009C10 
0C9C2C 
0C9C3C 
0090UC 
OC9050 
009C5C 
0090^0 
009080 
0C9C9C 
009  100 
0C91 IC 
00912C 
009130 
0C91*40 
0C915C 
0091GC 
0C9170 
0C9180 
009i90 


D 1 M L N S I C N  U(6,6) 
DC  I  J=l»6 
CO  ;   k=1,6 

1  l'(J,K)=0.C 
IF(  ;i-0)?,2,  ■ 

2  V=-1.C 


or 


IC.  U 


3  V=1.C 

1+  CT  =  CCSF(  THLrA) 
ST=SINF( TFETA) 
U(l  ,  1  )=CT 
b(  1  ,2  )=-V»ST 
U(2,  1 )=V«ST 
L(2,2)=CT 
U(3,3)=1.C 
U  ( U ,  u  )  =  1  .  C 
U(5,5)=lT 
U(5,6)=-*V«S1 
U(fc,5 )=V»ST 
U(6,6)=CT 
RETURN 
END 


SUPRCUriNF.  STIFCO  (  ThCT  A  ,RhC  ,  U  ) 
DIMENSICN  U(6,6) 
IF(RHC)  1  ,1,2 

1  V  =  -1.0 
GO  TC  3 

2  V=1.C 

3  CT=CCSr<THLTA) 
ST=SINF( THLTA) 
00    4     J=l ,6 

CO  4  K=l  ,6 
U  U(J,K)=0.C 
UM  ,1 )=CT 
U( 1 ,U)=ST*V 
U(2,2)=CT 
U(2,5)=ST«V 
U(3,J)=1.C 
U(l4  *  1  )=-ST»V 

u(u,u)=Lr 

U(5,2)=-ST«V 
U(5,5  )=CT 
U(6,6)=l  .0 
RETURN 
END 


SL'DPCUTIME    HANGlR     (CLX 
DIMENSION    U(fc,6) 
DO     1     J=l,6 
DC     1     K=1,6 


0092C0 

oc9:  10 

00922 C 
0092  50 
0C9c"0 
0C9.  ?C 
0C926C 
0091 t C 
0C9J  :0 
0C929C 
0095CC 

0  C  9  ;  "  C 

009320 
009250 

0C9S4C 
0C9350 
009i60 

C  C  9  I  "  o 
009330 

0C9.J?0 
0C9400 
0C9U1C 
0C9U2C 
OC9UiO 
0C9UUC 
0C945C 
009UeO 
0C9U70 
0C94°C 
009Uf)0 
0C95C0 
0C9510 
0C9-.20 
009  530 
009540 
009550 
0C956C 
00957C 
0C9530 
009590 
0096CC 
009olC 
009670 
0096 fO 
0C96UC 
009650 
0C966C 
009670 
0096E0 
0C969C 
0C9700 
009710 
009  72  C 
00973C 


1  U(J,K)=0.0  0C9/uC 

U(1,l)=l.r  C09750 

L(2,  2  J^1  .0  0C9760 

U(3,3)=1.C  CC9770 

U(»»,3)=CTZ  0C9f?C 

U(U,1*)  =  1.0  0C979C 

U(5,2)=-CLV  0C9i:00 

U(5f5)  =  l.C  0C9<:10 

U(6,1)=CLX  009t.C 

U(6,6)=1.0  009650 

RETC^N  0C9[;U0 

END  0C9£5C 

009260 
0  0  9  E  7  C 

SUBROUTINE  HANGlC  ( CT X , CL I , CTY  ,L )                                     009890 

PIMLNSICIS  U(6,6)  009F.9C 

DC  1  J=l,6  00990C 

DC  1  K^l  ,6  00991 C 

1  U(J,K)=0.G  009910 

U(1,1)=1.C  0C9910 

U(2,1)=C.TX  0099UC 

U(2,2)=i.C  0C9950 

U( 3, 5 1=1.0  0C99t)0 

IMU.U  )  =  1  '.0  009970 

U(5,U)=CTV        ,  0099=0 

U( 5 , J  >  =  1 .0  009990 

U(6,3)=-CL2  010CCO 

L(6,6)=1.0  010C1C 

RETURN  01CC20 

ENO  010C30 

O1OCU0 

0  ICC  5  0 

SUP3CUT1NE    PiONCH  (  SVf' ,  VV,  PH  ,  ,V  I  ,U  )                                                                                                OlCCoO 

ODlMcNSICNSVl3(6,22),  VV<  6t  6  )  ,U  (  6 , 6  )  t  R1  I  25,  251  »LL(25  )  ,l»f  (25  )  ,R2  (3,3  ) ,       010070 

1R(3,3),G1 ( 3,3) ,  C2< *,3 )  010C80 

v=i  010C9C 

DC  i  n=i  ,6  oioioc 

IF(SVB(N,M  ))  3,3,1  010110 

1  DO  _  L=l ,6  01012C 

2  VV(L,y)=VV(L,K)  •  0101 30 
M  =  V+1  0  1CUO 

3  CCNTINUc  0101'jO 
DO  H  L=l  ,3  •  0 10160 
Rl ( 1 ,L)=VV( 1 ,L)  0101 70 
R1(2,L)=VV(2,L)  01013C 
Rl (3,L»=VV( 3,L)  01C190 
R2( 1 ,L)=VV(M,L)  010200 
R2(2,L)=VV(5,L1  010210 

k    R2(i,L)=VV(6,L)  010220 

CALL  INVENT  (PI  ,?,0,LL,«M  010230 

CC  5  J=l,3  010240 

CO  5  K*1,3  OIO^C 

P(J,K)=C.C  010260 

DO  b    L=l,3  01C270 


J  R(J,K)=-  (J,K)+R.2(J,L)*R1  (L.K  ) 

SP=SINF(PHIJ  010230 

CP=CCSF(PHIJ  010^90 

G1(1,1)=CP  01C3C0 

G1M,2)  =  -SP  0  1  0 .-» 1  C 

G1(lfj)«C.  oici:c 

G1(2,1)=SP  0]G.-.iO 

G 1 C  2 ,  2  )  =  C  P  0  1  C  3  4  C 

Gl(2f^)=0.          •  0103-C 

Gl(i,1)=0.  Ol6i6C 

C1(5,2)=P.  010i70 

G1(3,J)-1.                                                .  010i?:C 

G2(1,1)»1.  010390 

G?(l,:)=0.  0  1  04  CO 

G2M,3)*o.  oiouio 

G2(2,U=0.  0104-0 

G2(2,2)=CP  010430 

G2(2,3)=-SP  010440 

G2l3,l)=0.  010450 

G2(5,?)=SP  01046C 

G2(i,i)=CP  010470 

00  6  J=l , Z  0  10480 

DO  C  K  = 1 ,3  010490 

R1(J,K)=0.0  010500 

HO  6  L=l,3  01  Ob  10 

6  Rl (J,K)=R1( J,K>+C2( J,L)«R(L,K)  010520 
OC  7  J  =  I  ,3  OlO'JiO 
OP  7  K=l , ?  01054C 
R(J,K)=C.O  010550 
DO  7  L=l ,3  010560 

7  R(J,K)=R(J,K)+RUJ,L).GHL,'<)  010570 
DO  G  J=l ,t  0105CC 
DO  8  K=l,6  010590 

3  U(J,K)=C.C  010600 

U(lfl)=l.C  010C10 

U(2,2)=1.0  010620 

U(3f3)=l.c  010630 

U(4,4)=1.0  010640 

U(5,5)=1.C  010650 

L'(6,6)  =  1.0  01066C 

DO  9  L=l,3  010670 

L(4,L)=R(1,L)  010680 

U(5f L)=R(2,L)  010690 

9  U(6,L)=R(3,L)  0107CO 

RETURN  0 107 1C 

END  010720 

010730 

0107UC 

SUBROUTINE  STflVEC(SV,N, BC)  01075C 

DIMtNSICN  RC(6,22)  0107cO 

DO  10  J=l,6  010770 

10  BC(J,NJ=C.  01078C 

IF(SV-1.)1, 1,2  010790 

1  BC(4,M=1.0  010800 

010P1C 


RC(3,N)=1 

0 

DC  (6,M  =  1 

0 

RE rURN 

if ( iv-r. > 

5 .  1 ,  U 

PC( 1 ,\)=1 

0 

R  C  ( 2  ,  N )  =  1 

0 

FCI3,N)=1 

0 

PLTL^N 

lF(SV-3. ) 

»S,6 

RC<3,N)=1 

c 

PClifN)=1 

0 

E  C  (  o ,  M  =  1 

0 

RETURN 

IF( SV-M .  )" 

'  r  7 »  e 

RC( 1,N)=1 

0 

hC(  .  ,.\)  =  1 

0 

BC  { 5 ,  N ) » 1 

0 

RCTUR'J 

BCf2,M  =  l 

0 

BC ( :  , N  )  =  1 

0 

RC(6,N)=1 

0 

BRANCC    (SV^.VV.PHIt1 
,U(6,6),VV(6 


M=1 

SR(  i  ,f I)=SVF  (  1 , y I 
SR(2,PI  »=SVH  ?,M 
SR(  5,1*1  >=SVt  (  U.,  M 

SR(»»,Mi)=sve(2»yi 

SRI  J,f  I  )=SVL'(  5,M 
SR(6,MI>=SVB{6,I"  I 
UV(  1,  1  )=VV(  1,1) 
UVMt2)*VV("l«3) 

UV(  1  ,3)=VV(  1 , k) 

uvt i ,ui=vv( l  ,2) 

UV(  1  ,5)=VV(  1,5) 
UV(  1  ,6)=VV(  1,6) 
UVU,1  )  =  VV(  3,  1) 
UV(2,2)=VV<3,2) 
UV(2,3)=VV( 3,U) 
UV(2,U)=VV( 2,2) 
UVU,5)  =  VV(3,5J 
UV(J,6)=VV( :,£ ) 

uvi  -:,  l  )  =  vv(u,  l ) 

UV(  3,2  )  =  VV(  U,  2  ) 
UV( 2,3)=VVCU,U) 
UV(3,4)=VV(U, 2) 
UV( i,5)=VV(U, 5) 
UV( S,61=VVU,6) 
UVU,1  )  =  VV(2,  1  ) 


(  3  ,  ?.  )  ,  I.  L  (  2  z  ) , 


I,UJ 

6  ) ,  V<  6 ,  *  > ,  R  1  (  25 , 2  5 ) 

w(25) ,SR( 6,22) ,LV(6. 


r^2(  3,  3),R(3,3), 


010P20 
0 1 C  e  3  C 
0  1  0  il  U  C 
0  10*50 
010060 

n  l  o  i:  7  c 

01CPP0 
0 10890 
010900 
010910 

01092C 
0  10930 
01C9U0 
0  10  9  5  0 
01C9O0 
01C970 
010930 
01099C 

0  1  1  C  C  0 
011C10 

01  1C2C 
0  11  C  3 0 
OllJuC 

onoso 

011C60 
0  1  1C70 
011050 

0  1 1090 

01  HOC 
011110 
011120 
Oil  130 
Oil  14C 
Oil  150 
Oil  160 
Oil  170 

omsc 

01  1  19C 
011200 
01121C 
0112^0 
011230 
01  12U0 
0112bC 
01  12cC 
011270 
0 1 12Q0 
01 1290 
0113C0 
01  lilC 
011320 
011.530 
0113U0 
011350 


uv(m  ,;  )=vvd  , 31 

UV( u , 3»=VV(?,U) 

uviM.ii  )  =  vv(  l, : ) 

l.V(4  ,:   )=WU,  51 

UV(t»,6)  =  VV{2,fr) 

UV(5, 1 )=VV(5, 1 ) 

UV(:>,2»  =  VV.(5,i) 

UV(oi  >)=VV( b  >  U  I 

uv(  ■_ », u  )-w( '  , : ) 

UVl 5  ,  5 )  =  VV( i, 5 ) 

UV(5f 6)=VV(5, 6) 

UV(c  , 1 )=VV(6,  1  ) 

UV(6,2)=VV(6,3) 

uvr c,7 )  =  vv( 6, u  j 

UV(6,U)=VV(6,2) 

UV(',,5)=VV(6,5) 

UV(6,6 )  =  VV< 6,6) 

DC    a    N*l,6 

IF(c,R(N,f'I  )  )  5,3 

1 

DO    I     L=1,6 

UV(L,n=UV(  L,M 

M*M»1 

CONTINUE 

no  u  k=t,3 

Kl  (  1  ,K)=UV1  1, K) 

Rl  U,K)  =  UV(2,K) 

Rl  (  S,K)=UV(  -:,k) 

R2(  1  ,,<l  =  UV|u,  K) 

R2 (2,K)=UV( 5, K) 

R2 (:,,K)=UV(6,K) 

CALL     INVERT     t  Rl 

3,0,  LL,KM 

DC    j    J=l,3 

cr    £    K=l  ,3 

R(J,K)=C.C 

DC    5    L  =  1  ,  5 

R{  J,K)=1( J,K) +R^ 

(J,L)»R1  (L 

SP=SINF(PHI  ) 

CP=LCSF(PHl J 

Gl  (1,1)=CP 

Gl(  1,2)=0. 

Gil  1,3)=-SP 

CI  12*1 1«0. 

6112*21*1. 

61(2,31=0. 

61(3,1)=SP 

G1  ii,l)=0. 

Gl  (io)=CP 

G2(  1  ,1 )=CP 

C2(  1  ,2  )  =  SP 

G2(  1,31=0. 

62(2,1)=-SP 

G2(-l,2)-CP 

62(2,31=0. 

62(3, 1)=0. 

U11570 
Oil  4  10 
01 1  i*7C 

01  luCO 
0  1  1  4  1  C 
0  1  1  u  i  0 
011«30 
01HHO 

0  1 1 4  5  C 
011460 
0111*70 

01  la~C 
01  U90 
011500 

0  11-310 
011520 
011530 

01  l-o  i«0 
01 1 55C 
01 1=60 
0  11370 
011530 

o  1 1 :  <i  c 

0116CC 
0  116  10 
0  1 162 C 

0  1  1  6  3  C 

01  16H0 

oii6:o 

0  1  1  6 1 C 
011c7C 
OlloPC 
0  1  1690 
011700 
0U71C 

0  11720 
011730 

01  17UC 
011750 
Ol17o0 
011770 
01 1780 
01179C 
Oil  POO 
011810 
011820 
0  11830 
0  1  1 8  a  0 
OllBbO 
011860 
011370 
0110P0 
011390 


S  (  J  ,  K  I  =  0 .  C 
L/C  7  L=1,3 
S(  J,K)=S( J, 

bo.  i    j=  i  ,< 


)  +P  1  ( J,L  )*C1  (L,K  ) 


DO 


=  1  ,> 


U( J,K>=O.C 
U  (  1 ,  1  )  =  1  . 0 

U  (  2  ,  1  )  =  ■j  (  1 ,  1  ) 

I  (  2  ,  2  )  =  1  .  0 

u < 2  , 5)  =  :> n  , : ) 

U(2,U)  =  sn,3) 
U  <  3 , 3  I  =  1 . 0 
U(4,U)=1.C 
UIS,U*S(2t  1) 

u  ( 5 ,  * )  =  5  ( : , ;. ) 

U  ( 5 , k I =S ( 2 , i ) 
U(5, 51  =  1.0 
U(6,l >=S(3,  II 
U  ( 6  f  2  I  =  S  ( 3 ,  7.  I 
U(6fU|=S(3,3l 
U { 6 . t   1=1.0 
RETURN 
END 


SUBROUTINE  STAVCC 
DIMENSION  BC(6,22> 
no  10  L=l,6 
10  BCfUN)=0.0 

I  F ( SV- 1.0)1.1,2 

1  B  C  ( 2 ,  N  >  =  1 . 0 
CC( 5,N1=1 .0 
RC(6,M  =  1.0 
RETUPN 

2  IF(SV-2.0)3,3,U 

3  BC{  1 , N  1  =  1.0 
BC(o,N)=1.0 
hC(U,N)=l .0 
RETURV 

U  IF(SV-3.C15,5»6 

5  BC(2»N)=1.0 
BC(U,N 1=1.0 
BC  (6,*J)  =  1  .0 
RETLRN 

6  IF( SV-U. 0)7,7,5 


0  1  l'JCO 

011010 
CI  19~C 
011010 
0  1  l')UO 

0  1  n  -j  o 

011960 

01  1970 
011-9S0 
01  1990 
012C0C 
012C10 
012020 
0  12C?0 
012CUC 
012CSC 
C12C60 
012C70 

oi2c  ;c 

012C9C 
012100 
012110 
012120 
012130 
C121HC 
0  121jO 
C1216C 
012170 
012180 
012190 
0122CC 
012210 
012220 
0122^0 
0122U0 
0  12250 
012^60 
0122  7C 
012280 
012290 
0123C0 
01251C 
012320 
012:5  50 
0123UC 
012350 
012360 
012370 
01233C 
012390 
0  12UCC 
012M0 
0  12*42C 
012U30 


fC(  I ,N)=1.C 
PCln  ,  M  -  1  .0 
PC ( 0  ,  H )  =  1 . C 


D!  V. 
SEA. 
0=1  .0 

nr-8(  k=i,n 

L  (  K )  =  K 

V(K)=K 

biga=a(k,k  > 

DP?0  I=K,N 

DC2C  J=K,K 

[  F  (  ARSF  (  f I G  A ) - AR  S F ( M  I ,  J  )  )  )  10,2,  ,20 

R I G  A  -  A 1 1  ,  J ) 

L!K)rl 

M  f  K  )  =  J 

CONl  INUL 

INT>  RCHANCE  RCfcS 

J=L(K) 

I  F  (  L  (  K  )  -  M  :  5  ,  3  5  , '  5 

C03C  1  =  1,  r; 

HCLL=-A(K, I ) 

A I K  ,  I  )  =  A  U ,  11 

A(  J,  I  )=H0LU 

INTERCHANGE    CCLL,vKf> 

I=M(K) 

I  F  (  f  (  K  )  -  K  )     ii  5  ,  4  .c  ,  ?  7 

DD4C    J=l  ,t-. 

HOLC=-A( J,K ) 

A(  J,K  )=,"  (J,  1  > 

A( J, [ )=hOL0 

DIVIDE   COLi:vN    ev    HVA    Pivcr 


1  = 


,  N 


CTN 
Lf'6 


I-K) 50,55,50 
,K)=A(I,K)/(-A(K,K)l 

IIN'U. 

ATRIX 


1  =  1  , 


<  ) 


,57 


IT ( J-K)    60,65,60 
A(I,J)=A<I,K)*A(K,J)  +  MI,J) 
CPNTIWJl 

01  VILE    ROW    bY    PIVTT 

DC75    J=l  ,N- 

IF  (J-K)  ,'0,7r. ,  7C 

A  (  K  ,  J  )  =  A  (  K  ,  J  )  /  A  (  K  ,  K  1 


012470 
0  1  £  K  '  0 

0  )  2  4  0  C 

U12.O0 
0 12.10 

o  1 2  -. ;  c 

0  1  2  'J  5  C 
0125*0 
012550 
0  1 2  -j  A  0 
0125/C 
0  12-j."'C 
0  12'jOO 
0126C0 
Cl2ol0 
0126  20 
012630 
012640 
012650 
0  12060 
0  1  2c 70 
01268C 
012000 
012700 
01271C 
0127"C 
0  12730 
012740 
0  1270 
0  12760 
012770 
0  12730 
012700 

oi2rcc 

012810 
012320 
01 2830 

o  1 2  ;■  i*  o 
012350 

012o60 
012370 

oucc 

012890 
012900 
0 1291 C 
01292C 
012930 
012940 
0129jC 
012960 
012970 


CI  Ml  IMJi 

com i\Ul 

i   p-i  cue  t  nr  m  i vr  t s 

D=0«A(K, 

K  ) 

R[   PLACF 

PIVCT  HY  KCCIPRCCM. 

4  (  K  ,  K  )  =  1 

.(  /  M  K,  K  ) 

C.TN1  1*  [):. 

FINAL  "!■: 

:V»  ANl  CCLI  I*-?  INTEf CI 

K  =  N 

K= { K-  1J 

t  r  (  h  )  !  . 

0,  150,  IC- 

I  r  { i  -  k  ) 

i?o ,  i :o,  it s 

DC1  IC  J  = 

1  iN 

HPLi =A(. 

l»K  ) 

M  J  ,  K ) =- 

A  (J,  I  ) 

A{  J,  I  )=l 

IC  LU 

J=M(K) 

[ft J-K) 

100,  ICC,  125 

DO  130  1= 

1  ,N 

HCLL=A(K 

,  1  ) 

MK,I)=- 

M  J,  1  ) 

A( JfI l-H 

■OLD 

on  re  lc 

10 

RETURN 

SUBRCUT  INE  TINNAT  (U,UU,  I,V) 
PIMCNSIOf  U(t  ,6  )  ,L'U(g, 6  )  ,V( £  ,6  ) 
I  F  (  I  -  U  t   1  ,  6  ,  6 1 
62  l/C  ?C  J=1,A 
DC  ,'0  K=l  ,6 

70  V( J,K)=C( J,K) 
RETURN 

61  CO  71  J=l,r. 
DC  ,'1  K=l,6 
V  (  J  ,  K  I  =  i  .  C 
CO  71  L= 1 ,6 

71  V(J,K)=V(J,K)+U(J,L >*UU(L,K) 
RETURN 

bNL) 


SUBROUTINE  F1NBRA  (  U,  VV,  MUV< ,  V  ) 
DIMENSION  U(t,6),VVl6,6),VU,6) 
IF INNHR-  I  11,1  ,3 


i  dc  :  J= 
DC  I   K=i,e 

2  V(J,K)=UU, 
RETURN 

3  DC  u  J=l  ,t 
DO  h    K=l  ,6 

V  (  J  ,  K  )  =  C  .  C 
DO  i4  L  =  l  ,t 


oi3c: 

;o 

CI  50 

!  0 

o  1 3  c : 

:c 

0  131 

50 

0  1  30' 

01  3C: 

lc 

0  13O 

-0 

0 1  3 1 C  0 

3  13  110 

o  i  i  i  :>  c 

013130 

0  131UC 
C1315C 

0 1  3  1  6  C 
D13  i  'C 
013130 
013100 
013200 
0  1  3 1  i  C 

013220 

013_i0 
0132'iO 
013,  5G 
013<_60 
013^70 
C13^°C 
0  13290 
013300 
013310 
01 33? C 
013330 
0135UO 
013350 


0  1  I 


c3  0 


0  13390 

01  3U00 
013U10 
O13U20 

0  1  3*4  30 
013UUO 
013U5C 

01  3^oC 
013^70 
C13VC 
0  13U9C 
013500 
0  1  3  5  1  C 


V(J,K)=V( J,K)+U(J,L)«VV(L,K) 
RETLRN 


SUB30UTIME  CELPA(LU#SVI,SVE.V) 

DIMENSION  UU(6,6)',VV(6,6),V<6,6),SVI<6,1),SVEU,1) 


DC  9C  N=i,6 

IF ( S VI (N,l)) 1,90,1 
1  K=M+1 

DC  91  J=l  ,6 

91  VV(J,M)=UU(J,N) 
90  CON  II  NUT 

L=C 

DC  9?  N=1,6 
IF(SVE(N, 1 ) 192,2,92 
2  L=L+1 

DP  9?  K=l ,3 
93  V(L,K)=VV(N,K) 

92  COM  1NUE 
RETURN 
END 

END 


013S2C 
013530 
0135UC 
C13550 
O13'j60 
013570 

0  135  -C 
013590 

01  3600 

o 1 3e  ic 

01362C 
0  13630 
0136U0 
0136SC 
01366C 
013670 
013680 
013690 
0137C0 
013710 
013720 
013730 
0137^40 
013750 


APPENDIX  e 
ACCURACY  OF  METHOD 


E.l   General 

Appendix  E  contains  the  vibration  analyses  of  the 
several  typical  piping  configurations  to  which  reference 
is  made  in  Chapter  3<»   The  analyses  are  grouped  according 
to  the  purpose  for  which  they  were  performed.   The  firs; 
group  was  performed  to  establish  the  accuracy  of  the  method 
embodied  in  the  program  VIPIPE,  the  second  and  third  to 
establish  its  integrity. 
E.2   Method  Accuracy  Analyses 
E.21   Systems  Analyzed  (straight  sections) 


System 

Parameters 

Systems 

1   *<    i 

\          2' 

\       3* 

end  conditions 

\                      \ 
fixed-fixed 

^fixed-free 

f ixed-oropped 

length  (in.) 

200 

200 

200 

diameter  (in.) 

2.0 

2.0 

2.0 

wall  thickness  (in.) 

.125 

.125 

.125 

density  (lb./ou.ft.) 

1172 

1172 

1172 

shear  modulus  (psi) 

30  •  106 

30  '  106 

30  e106 

elastic  modulus  (psi) 

12°  106 

12-  106 

12=  106 

E-l 


o 2 2   Comparison  of  Natural  Frequencies 
a0   System  1 


Mode 
number 

frequency  (rado/secJ 

Difference  {%) 

VIPIPE 

Comparison  Values 

1 

48«53300867 

48,533126 

0 00024 

2 

133°78346278 

133*78345 

oOOOOO 

3 

262.2684309? 

262,26886 

,00016 

4 

433*54478122 

433 o54393 

0OOOI9 

5 

64? 062290622 

647 o63991 

O0026 

b0   System  2 


Mode 
number 

frequency  ( rad o/sec) 

Difference  {%) 

VIPIPE 

Comparison  Values 

1 

7 o62356485 

7o6270709 

c 04595 

2 

47 o79657598 

47o79776l 

o002^-7 

3 

133o8363688 

133o83635 

oOOOOl 

4 

262o26528l07 

262026544 

c00006 

5 

^33 0 5^435149 

433 0 54415 

0 00004 

System  3 


Mode 
number 

frequency  (rad0/sec0) 

Difference  {%) 

VIPIPE 

Comparison  Values 

1 

33 o44525745 

33o445600 

0 00102 

2 

108038602303 

108c38599 

0 00003 

3 

226ol3877l5 

226ol3868 

0 00000 

4 

386071069206 

386o71057 

o00003 

5 

590,10235871 

590,10164 

0 00012 

Sources  of  comparison  frequencies  for  problems  1 
through  3  are  given  in  Ec23& 

E»2 


E.23   Source  of  Comparison  Frequencies  for  oysters  1  through  3 

The  governing  fourth  order  equation  of  a  homogeneous 

single  component  straight  section  system 

viv  -  na>2 
y  EJ~* 

was   solved  to   give 

y   =  a-cos^x  +  b«coshocx  +  c«sinp(x  +  d«sinh<*;x 

Upon  applying  boundary  conditions  for  problem  1 

(y=yf  =0  at  x=0,  and  x=L)  and  simplifying  it  was  found 

that  the  eigenvalues  of  the  system  must  satisfy 

cos oc  co shoe  =1  (E.23.1) 

where 


V  EJ 

Upon  applying  the  boundary  conditions  of  problem  2 
(y=y'=0  at  x=0,y"=y",=o  at  x=L)  it  was  found  that  the 
eigenvalues  of  the  system  must  satisfy 

cosoCooshcc=-l  (Eo23-2) 

Upon  applying  the  boundary  conditions  of  problem  3 
(y=y'=0  at  x=0  ,y=y*=0  at  x=L)  it  was  found  that  the 
eigenvalues  of  the  system  must  satisfy. 

sinoccoshoC  -cosocsinhcc=0  (E823~3) 

The  roots  of  equations  E«23-l,  E023-2  and  E.23-3,  and 
from  them  the  first  five  natural  frequencies  of  the  systems, 
were  obtained  through  the  use  of  a  digital  computer • 
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E92^   System  Analyzed  (section  of  constant  curvature} 
System  parameters s 


radius  of  curvature  (in8) 

50 

included  angle  of  arc 

2?0° 

diameter  (in*) 

2 

wall  thickness  (in0) 

0125 

elastic  modulus 

30  °  10°"  psi 

shear  modulus 

12  •  106  psi 

density  (lbs  0/cu«-ft «, ) 

556 

Results: 
in-plane  vibrations* 


end  conditions  %      fixed-fixed 


Mode 
number 

frequency  (rado/seco) 

Difference  {%) 

VIPIPE 

Comparison  Values 

1 

73o329730?2 

73*6 

o23 

out-of=pIane  vibrations 


Mode 
number 

frequency  (radu/sec0) 

Difference  {%) 

VIPIPE 

Comparison  Values1 

1 

38o87735138 

39o8 

2o3 

E°3   Method  Integrity  Analyses 

The  first  group  of  analyses  considered  here  constitutes 
a  mirror  image  check  on  system  integrity 0   Four  configura- 
tions were  considered,  and  several  mathematical  models 
were  used0   The  second  group  consists  of  two  analyses  of 


I 
This  value  was  obtained  usin^  a  curve  in  reference  2 
which  'was  constructed  through  the  use  of  a  modified  Rayleigh 
method  •   In  the  words  of  the  reference  5  "  0 «.  0 laboratory  tests 
have  shown  the  curves  to  be  substantially  correct 0" 

Comparison  values  obtained  by  employing  the  equation 
CO  =F(qQ  2  fET}   after  obtaining  F   from  curves  In  ref0  20 

E-4 


the  same  single  branch  system  v/herein  first  one  member  then 
another  is  considered  to  be  the  brancho 
E.31  Mirror  Image  Analyses 
a0   System  1 


system 
mirror  image 


System  parameters:   See  system  2  E„21 

Mathematical  model:   Distributed  mass  with  the 

effects  of  shear  deflection  and  rotational  inertia  considered. 

Results: 


Mode 
number 

frequency  (rado/sec.) 

Difference  (>2) 

System 

Mirror  Image 

1 
2 
3 
4 
5 

7.622895556 

47.76681871 

133»636l4681 

26I.54431793 

431.64844506 

7o62289556 

47.76681872 

133»636l4691 

261 o54431024 

431 o64917860 

,00000 
.00000 
o00000 
,00000 
,00016 

b9   System  2 


system 
mirror  image 


System  parameters:   See  system  3»  E.21 
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Mathematical  model:   Distributed  mass  with  the 
effects  of  shear  deflection  and  rotational  inertia  neglected. 

Results: 


Mode 
number 

frequency  (rad,/sec) 

Difference  {%) 

System 

Mirror  Image 

1 

33.^5257^5 

33.^5257*5 

.00000 

2 

108o38602303 

108o38602303 

.00000 

3 

226ol3877715 

226,13877715 

.00000 

4 

386,71069206 

386.71059^81 

.00002 

5 

590.10235871 

590.08973387 

.00211 

c.   System  3 

i 1 


system  (component 
numbers  shown) 


indicates  a  pipe 
hanger 


►mirror  images 


System  parameters: 

Component 

Parameter 

1 

2 

3 

k 

5 

6 

7 

diameter  (in.) 

2.0 

2,0 

2,0 

H 

2.0 

2,0 

2,0 

wall  thickness  (in.) 

.125 

.125 

.125 

N 

0125 

,125 

.125 

length  (in.) 

130. 

1,^0 

6,0 

G 

60.21 

1,50 

6,28 

radius  of  curv{in.) 

CO 

ij-.O 

00 

E 

30. 

cx=> 

fc.O 

incl.  angle  of  arc 

0 

20° 

0 

R 

115.° 

0 

90° 

E-6 


Intensive  properties  of  all  components  except 
hangers; 

density:   556  lbs./cu.-ft. 
shear  modulus:   12  •  10°  psi 
elastic  modulus:   30-  106  psi 


Hanger  spring  rates: 

CLX  =  1000  lbs. /in. 

CLI  =  CLZ  =  1000  lbs. /in. 

CTX  =  CTY  =  CTZ  =  1000  in.-lb./rad. 

End  conditions: 

Both  ends  of  system  fixed. 

Mathematical  model:   Distributed  mass  with  the 
effects  of  shear  deflection  and  rotational  inertia  considered. 

Results: 

in-plane  vibrations 


Mode 

number 

frequency  (rad./sec.) 

Difference  ($) 

System 

Mirror  Image  a. 

1 
2 
3 
k 
5 

116. 91660080 
305.21721131 
567.56373058 
8^5.99^32^21 

116.9166008^ 
305.21720649 
567o56397808 
846.1762^950 
1135.09717500 

.00000 
.00000 

.00005 

.02150 
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out -of -Diane  vibrations 


node 

number 

f req uency  ( rad  .  /sec  . ) 

Difference  {%) 

System 

Mirror  Image  a0 

1 
2 
3 
i* 

5 
6 
7 

114.21818694 
2^3-92831220 

321c'4089H23 
591.40446121 

872.23102583 
1138o01776773 
1571«78173050 

114=21818694 
243.92831317 
321.40891770 

591 0 40449774 

872 0 24253187 

1138ol96l5015 

,00000 
,00000 
,00000 
,00000 

.00131 

c 01567 

in-plane  vibrations 


Mode 
number 

frequency  ( rad ,/sec . ] 

Difference  {%) 

System 

Mirror  Image  b. 

1 

2 
3 

4 

Il6o9l660080 
305.21721131 
567o56373058 
845 • 99432421 

Il6o9l660080 
305o21?21131 
567056373058 
845 o9943421 

*00000 
oOOOOO 
000000 
,00000 

out-of -plane  vibrations 


Mode 
number 

frequency  (rad./seco) 

Difference  {/<>) 

System 

Mirror  Image  bo 

1 

114o21818694 

114o21818694 

.00000 

2 

243e92831220 

243c9283i220 

oOOOOO 

3 

321o4089H23 

321,40891123 

.00000 

4 

591.40446121 

591 0 404461 2 1 

.00000 

5 

872.23102583 

872.23102583 

.00000 

6 

1138c01776773 

1138c01776773 

oOOOOO 

7 

1571.73173050 

1571o78173050 

oOOOOO 

£=8 


system  (component 
numbers  shown) 


J&  indicates  pipe 
hanger 


airror  image 


System  parameters: 


Component 

Parameter 

1 

2 

3 

k 

5 

6 

7 

diameter  (in.) 

wall  thickness  (in.) 

length  (in.) 

radius  of  curv.(in.) 

incl.  angle  of  arc 

2.0 
.125 

^00. 
0° 

2.0 
.125 

39.27 

50. 

4-5° 

2.0 
.125 

100. 
0° 

2.0 

.125 
39.27 
50. 
45° 

H 
A 
N 
G 

2.0 

.125 

100. 
0° 

2,0 

.125 
100. 

CXD 

0° 

Component 

Parameter 

8 

9 

10 

11 

12 

13 

diameter  (in.) 

wall  thickness  (in.) 

length  (in.) 

radius  of  curv.(in.) 

incl.  angle  of  arc 

2.0 

.125 

39.27 

50. 

45° 

H 

A 
N 
G 
E 
R 

2.0 

.125 
39.27 
50. 
k5° 

2.0 

.125 

100. 

0° 

2.0 

.125 

39.27 

50. 

V>° 

2.0 

.125 

100. 
0° 
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hangers : 


Intensive  properties  of  all  components  except 


density?  k$Q   lbs0/cu0-fto 
shear  modulus:   12°  10  psi 
elastic  modulus  30  -  10°  psi 

Hanger  spring  rates: 

CLX  =  CLY  =  CLZ  =  1000,  lbso/in0 
CTX  =  CTY  =  CTZ  =  1000o  ine~lb0/rad0 

End  conditions:   all  aiain  member  and  branch  ends 


fixedo 


Mathematical  models   Distributed  mass  v;ith  the 
effects  of  shear  deflection  and  rotational  inertia  considered. 

Results : 

in-plane  vibrations 


iviode 
number 

frequency  ( rad/s ec „ ) 

Difference  (/&) 

System 

Mirror  Image 

1 
2 
3 
4 

18,55192182 
51*25211856 

I8o55l89208 
51.25212736 

90,9^757587 
107ol9459215 

,00016 
,00001 

E-10 


out-of-plane  vibrations 


node 
number 

frequency  ( rad ./sec.) 

Difference^) 

System 

Mirror  image 

1 

14.31754957 

I4c31754879 

.00000 

2 

3^,45495449 

3^^5^9/4654 

e00G02 

3 

41.07914530 

4lc07913499 

.00002 

4 

73 0 06452245 

73.06429798 

.00030 

5 

115.86572509 

115o86505834 

.00057 

£.32   Single  Branch  System 


"lain  Member 


Case  1  (component 
numbers  shown; 
(p  ~180° 


Main  Member 


System  Parameters: 


Case  2      (component 
numbers   shown) 
0  =  180° 


Component 

Parameter 

1 

2 

3 

4 

diameter  (in.) 

wall  thickness  (in.) 

length  (in.) 

radius  of  curv.(in.) 

incl.  ancle  of  arc- 

2.0 

.125 

100. 

0° 

2.0 

ol25 

100, 

0° 

2.0 
.125 
,01 
.01 

45° 

2.0 
=  125 

100. 

0° 

£-11 


Intensive  properties  of  components: 

density:  55&   lbse/cu0~ft0 
shear  modulus:   12-  10  psi 
elastic  modulus:   30"  106  psi 

End  conditions:   Both  ends  of  main  member,  and 
end  of  branch,  fixed . 

Mathematical  model:   Lumped  mass  with  the  effects 
of  shear  deflection  and  rotational  inertia  neglected. 

Results: 

out-of-plane  vibrations 


Mode 
number 

frequency  (rado/seco) 

Difference  {%) 

System  1 

System  2 

1 
2 
3 
if 
5 
6 

^9»6^60^938 
135 05268^605 
269o83567590 
530o720123l6 

675o5l662199 
901.52786818 

4-9  O6^60^96i 
135o52685066 
269o83569992 
530o757l8901 
675o570^008 
901 076873827 

oOOOOO 
oOOOOO 
oOOOOO 
c 00698 
o00796 
0 02671 
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In-plane 

vibrations : 

Mode 
number 

frequency  (rad./sec.) 

Difference  {%) 

System  1 

System  2 

1 

132,89829688 

132,39836401 

.00005 

2 

190.61015386 

190.61009901 

.00002 

3 

^30.03806330 

430,03948554 

.00033 

4 

512.55370411 

512,54463160 

.001?? 

5 

893o66594383 

893.10993095 

.06221 

6 

962,07559083 

963800709091 

,09682 

E-13 


APPENDIX  F 
TRANSFER  MATRIX  oOUrlCE  REFERENCE 


F.l 


General 


Fourteen  of  the  transfer  matrices  used  in  VIPIPE 
were  drawn  from  the  catalog  of  transfer  matrices,  pages 
376-MO,  in  Pestel's  "Matrix  Methods  in  Elastomechanlcs" 
DG  .   Some  of  the  transfer  matrices  constructed  in  VIPIPE 
are  a  composite  of  two  separate  transfer  matrices  f  ro  a  [1] 
These  composite  matrices  are  constructed  as  shown  in 
Figs.  F.l-1,  and  F.l-2.   The  former  shows  the  format  of 
case  1  composite  transfer  matrices,  the  latter  those  of  Case  2. 
In  both 


'1,1 


>2,1 


lj.2 


2;2 


null 
submatrice 


Fig.  F.l-1   In- plane  composite 
transfer  matrix  format 


cea 


Fig,  F.l-2  Out-of -plane 
composite  transfer  natrix 
format 


cases  the  a  matrix  is  of  fourth  order,  the  B  natrix  of  second 
order.   In  Case  1  composite  transfer  matrices  the  elements 
of  the  B  matrix  occupy  the  positions  shown  (Fig.  F.l-1). 
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In  the  following  source  listing,  when  two  matrices  are 
indicated  for  one  subroutine,  the  first  is  the  A  matrix 0 
(All  component  transfer  matrices  are  constructed  in  sub- 
routines 0  ) 

F o 2   In-Plane  Vibration  Transfer  Matrices 
Transfer  matrix 


straight  section  with  distri- 
buted mass,  shear  deflection 
and  rotational  inertia 

straight  section,  elastic 
massless  field  with  shear 
deflection 

curved  section,  elastic  mass- 
less  field  without  shear 
deflection 

point  mass  with  rotational 
inertia 

rigid  straight  section  with 
rotational  inertia,  distributed 
mass 

stiff  corner  (elbow),  massless 
field 


Subroutine 

Source 

DI3T14 

C-3a 
C-la 

SFIELD 

C-3c 
C-lo 

CFIELD 


POINT 


RIGID 


STIFCO 


* P . 3  Qut-of-Plane  Vibration  Transfer  Matrices 

Transfer  matrix  Subroutine 

straight  section  with  dlstri-    DISTrtO 
buted  mass,  shear  deflection 
and  rotational  inertia 

straight  section,  elastic  mass-  SFIEIO 
less  field  with  shear  deflection 

curved  section,  elastic  mass-    CFIELQ 
less  field  with  shear  deflection 


point  mass  with  rotational 
Inertia 

rigid  straight  section  with 
rotational  inertia,  distributed 

mass 

stiff  corner  (elbow),  massless 
field  F-2 


C-?c 


C-3f 
C-lb 

C~3<i 
C-lb 


C-?f 


Source 

C-3a 
C~2a 


C-3c 
C-2c 

C-8c 


POINO 

c-3f 

C-2b 

RIGIO 

C-3d 
C-2b 

STIFOO 

C-8f 

APPENDIX  G 
SAMPLE  PROBLEM 

G»l   General  Remarks 

Appendix  G.  is  intended  to  supplement  appendix  B, 
Included  in  appendix  G.  is  an  illustration  of  the  data 
deck  for  a  typical  multiple  branch  system  with  hangers, 
and  an  example  of  the  three  output  formats  available  in 
VIPIPE. 
G.2  The  Data  Deck  (f  Ug,  Q.2rZ) 

The  data  deck  as  shown  calls  for  determination  of  both 
Case  1,  and  Case  2  natural  frequencies  (IoP=0).   If  only 
Case  1  frequencies  were  sought,  card  ^  would  be  removed 
from  the  deck,  and  IOP1  (circled  field  on  the  first  data 
card)  would  be  set  equal  to  1»   Similarly,  if  only  Case  2 
frequencies  were  sought,  cards  2  and  3  would  be  removed, 
and  IOP  set  equal  to  20 

In  the  sample  problem  configuration  all  boundary  con- 
ditions may  be  specified  by  code  (B86l),  but  in  order  to 
illustrate  the  method  of  programming  boundary  conditions 
for  which  no  code  is  provided,  the  Case  1  boundary  condition 
state  vector  for  branch  1  is  read-in  directly „   Note  that 
the  second  field  of  card  2  is  blank  (blanks  are  read  as 
zeros  by  the  computer)  and  that  card  3  contains  the  particular 


dee  a open,  B.,  section  B*71« 


G-l 


BiMiiCK 


4/^C 


Co.,, 


Branch         2. 


!pr-»    ^-f 


Flo.     G-.2-1         Multiple    Branch    Svstevn 


460. 

30.00E  06  12.00E 

06 

87 

".  1 25      1 00  7' 

7 125  -50.     "45. 

i     2« 

.125      100. 

a. 

.125  50.       45. 

2000.   1000. 

8. 

.125  -50.       45. 

2000.   1000.   1000.   1000.   1000.   1000. 


2.   .125  50.        457" 

~~ T.      ^125  100. 
~27~  :if5~5Q~.    ~  45. 

~Ti      TT25  400." 

T2fr"3E~~ — 

1    l'  *•  *• 

1.1.1. 


I 


•       II       II  I 

|o|ooooo|oo|ooo||lo(iiiuo|ooouaoooi 


'222222221222 
I3|3333  33  3|3  3 


2222222 
333|333 

313  313  3 

)3|3 

5555555 

5555551 

5  5  5  5 

666(666 

6666666 

6  E  b  E 

8S8IMS 
9999999 

9  9  9  9  9  9  9 

9999 

5555555555555 
6666666666666 


19  9  9  9  9  9  9  9  9  9  9 


lllllllllii^tll 
222222 


155  5555 
i  6  G b 6 66 


3  .  5  I  o  B  8  ! 
9999  j999 


333333333333 


555i5bob5555 
66b  6  b  -,  C  ?  f  666 


I  ! 


.2-2    ......  ■  |  I.  Datn'l  A 

C--2 


boundary  condition  state  vector e 

The  data  deck  used  in  Fig.  G02-2  pertains  to  system  h 
of  section  E.31,  appendix  E0   The  system  is  illustrated  in 
Fig.  G.2-1. 
G.3   Output  Options 

Examples  of  the  three  output  formats  available  in 
VIPIPE,  through  specification  of  the  value  of  the  parameter 
GijxiBC  (B.7D,  are  included  on  pages  G-4,  and  G-5o   The 
system  analyzed  was  system  3  of  section  E,31»  appendix  E. 
and  is  illustrated  there.   The  problem  numbers  refer  to 
the  order  in  which  the  data  decks  are  "stacked"  (B.8). 
The  heading  ...  "PROGRAM  VIPIPE"  . ..  etc.  is  printed  once 
per  run,  e.g*,  if  ty  data  decks  were  "stacked"  behind  the 
program  deck,  and  each  asked  for  Case  1  and  Case  2  analyses 
the  heading  would  precede  the  print  out  of  the  results  of 
the  Case  1  analysis  of  system  1  only0 
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PROGRAM  VIPIPE 

METHOD  OF  TRANSFER  MATrTcES 


MOCAL  FRECLENCIES  D 


PROBLEM.  TERMINATED  CLE  10  SIGNIFICANT  FIGURE  LIMITATION  OF 

A  OISTRIBUTEC  MASS  APPROACH  WAS  EMPLOYED  AND 
THl  EFFtCTS  OF  SHEAR  DEFLECTION  ANO  ROTARY  INERTIA  hERE  CONSIDERED. 


PIR  SECCN 
LENCY 
6.9166008M 

5.J1T2C89: 

7.94W9170 

6.1U0U5651 

CCMPUTER. 

Fig.  £3-1        GDABC  =0. 


FRCtRAM  WF1PE  C  E.  FHM   NHA-1  FEBRUARY  Hit  . 

MODAL  FREQUENCIES  3F  A  PLANAR  PIPING  SYSTEM  ARE  DETERMINED  BY  AN  ITERATIVE  PROCEDURE  USING  THE   ' 
METHOO  PT  TRANSFER  MATRICES. 


PROBLEM 

1 

IN  PL/WE  MODE  FRE8UENC 

ES  (RAOIANS  PER  SECCNDI 

NOOE 
1 
2 

3 

5 

FREQUENCY 

lit. 91640080 
.  Iff, 21721073 
5*7. 5*385876 
81A.05107199 
1128.3763U979 

PROBLEM  TERMINATED  DUE  TO  SIGNIFICANT  FIGURE  LI 

I  TAT  ION  OF  CCMPUTER. 

THE  EFFECTS  OF  SHEAR  DEFLECTION  AND  ROTARY  INERTIA 

E.imfMH. 

HALL  TMCKNESS      LENGTH 


*q 


mm    anclI^fIse 

•49  20:ll 


!  ■ 

IODULUS 
i|Jj 


naggt&s 


Ml.-   :\l 


JUNDARY  CONDITIONS 


Fig.  G.  3-2 


GDABC  =1. 
6-4 


i  PER  SECCNCI 


PROBLEM  TERI 
THE  EFFECTS  OF  SHEAI 


DEFLECTION  ANC  n 


FIGURE  LIMITATION  OF 
)ACH  WAS  EMPLCYtC  AND 
:rtia  WERE  (  ' 


SECTION  LENGTH  ANO  SLBSECTION  CAI 


ITERATION  NUMBER 


Mfgf 

r. --.■  :  •  9 


...■■TIL'  It 

I 

7679791  " 

| 


I 


.  ,,■■■  ■ 


- 

hi  -    -    r 


:8t 


556. OC 

556. OC 

illi 


.30E»OB 
.30E+08 
.3CE«08 
.3OE»08 


Fto.  G.3-3 
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GDAB0  3. 


